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C*i Abstract 

O _ , _ 

We investigate rare semileptonic B — > K*£ + £ decays, providing a 

comprehensive treatment of theoretical uncertainties in the low-g 2 region 
as needed for interpreting current and future LHCb and 5-factory data in 
i terms of the new physics search. We go beyond the usual focus on form- 

factor uncertainties, paying proper attention to non-factorizable terms. 
A central point is the systematic exploitation of the V — A structure of 



Oh 



Qh SM weak interactions, which leads to the suppression of two helicity am- 

plitudes and some of the angular coefficients. We review how this works 
'— 1 at the level of (helicity) form factors, and show that the hierarchies extend 

^ to non-factorizable terms. For virtual charm effects, we give an argument 

for it in terms of light-cone QCD sum rules that continues to hold at the 
level of "long-distance" Aq CD /?ti 2 power corrections, reducing an impor- 
ts) tant source of theoretical uncertainty in any B, B s — > Vl + l~ (or B — > V7) 
decay. The contributions of the remaining hadronic weak Hamiltonian re- 
CN spect a similar hierarchy. We employ a resonance model to preclude (in 

the B —7- K* case) large long-distance corrections to this, 
^—i A phenomenological part pays particular attention to the region of 

lowest dilepton mass, 4m| < q 2 < 2 GeV 2 . Two observables remain theo- 
: '"j retically clean, implying a (theoretical) sensitivity to the real (imaginary) 

part of the "right-handed" Wilson coefficient C' 7 to 10% (1%) of Cf M , both 
in the muonic and the electronic mode. We also show that there are two 
near-exact relations between angular coefficients, even in the presence of 
new physics and when lepton masses are not neglected. 
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1 Introduction 



Rare B decays are CKM and/or loop-suppressed in the Standard Model (SM), 
which accounts for their small branching fractions and makes them sensitive to 
small contributions from possible new degrees of freedom beyond the SM (BSM). 
Unfortunately, they generally suffer from non-perturbative strong-interaction ef- 
fects, the perennial bugbear of flavour (and collider) physics and difficult to con- 
trol theoretically. An example where such uncertainties are small is provided by 
the decay B s — > for which LHCb has just published first evidence [l], 

BF(B S -)■ ^V)e xp = (3.2i^| stat l°!| syst ) x 1(T 9 [3.5a significance], (1) 

to be contrasted with the SM theory prediction, 

BF(B S fx + fi~) S M = (3.54 ± 0.30) x 10~ 9 , (2) 

which incorporates estimates of soft-photon emissions and takes account of the 
finite lifetime difference of the B s [2j[3]. This result strongly constrains new- 
physics contributions through scalar x scalar semileptonic operators, such as may 
occur in the MSSM at large tan/3. It also constrains physics entering through a 
modified Zbs vertex, although less strongly, and in view of the statistics foreseen 
at LHCb, this will remain so for the foreseeable future. Beyond this single, 
exceptionally clean example, one generally has to deal with amplitude ratios 
including strong rescattering phases, which are currently not tractable with a 
first-principles method like lattice QCD. One faces a trade-off between theoretical 
cleanness (as in leptonic decays) and number of modes, which increases with the 
number of hadrons in the final state, as does the theoretical complexity. 

Rare semileptonic and radiative decays such as B — > K*fi + fi~ and B — > K*^ 
provide, in a sense, the best of both worlds. On the one hand, they have a 
rich kinematic and helicity structure. Experimental results on B — > K*fi + fi~ 
include measurements of the forward-backward asymmetry and the longitudinal 
polarisation fractions by the -B-factories, CDF, and LHCb [4T[7], of the angu- 



lar observables and A im by CDF |5j 8] and related angular coefficients by 
LHCb j6], and measurements of non-angular observables (e.g. (9||TT]). Results on 
the full angular distribution are anticipated from LHCb. B s — > has also 



been observed 12 13 . On the theory side, there have been a large number of 



conceptual and phenomenological studies (e.g. 14-68]). On the other hand, in- 
volving only one hadron in the initial and final states, these modes approximately 
factorize naively, into products of nonperturbative form factors and Wilson coeffi- 
cients given in terms of short-distance SM and BSM vertices. Naive factorization 
is, however, violated by the hadronic weak hamiltonian. Schematically, 



A(B^V£-£ + ) = J2 C i( ri+ \ lT i l \°)( V ^ r 'i b \ B ) 



+^<^|i^|0>i ? .r.<V|T{jSi(ar)«S d (0)}|B> I (3) 
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where C» are semileptonic Wilson coefficients, is the hadronic part of the 
electromagnetic current, and is the hadronic weak AF = 1 hamiltonian. 

This equation is correct up to higher orders in the QED coupling a em . The 
form factors (V|sr^6|JB) are, in principle, accessible on the lattice. For small 



dilepton invariant mass, they also lend themselves to light-cone sum rules 69 



and satisfy certain constraints from heavy-quark symmetry/large-energy relations 



70 . Making use of these relations, a number of observables can be shown to have 



reduced sensitivity to hadronic (form factor) uncertainties (see eg 15 17 32 57]). 
However, this cleanness is compromised by the non-factorizable terms on the 
second line of ([3]), which pose the most challenging theoretical problem. 

A breakthrough has been a comprehensive treatment of these effects within 
QCD factorization (QCDF) based on the heavy-quark limit 23 24 , valid for large 
hadronic recoil and in particular in the "low-g 2 " region of dilepton invariant mass 
below the charmonium resonances, and for the radiative decay. At leading power 
in an expansion in A/mb, where A = Aqcd is the dynamical QCD scale, the terms 
on the second line in ^ are expressed in terms of form factors, decay constants, 
and convolutions of light-cone distribution amplitudes of the initial and final 
state mesons with perturbative hard-scattering kernels. The same framework also 
allows one to systematically compute perturbative corrections to the form factor 



relations 71 . QCD factorization breaks at subleading powers, due to end-point 
divergences. The power corrections that do not factorize have to be modeled or 
constrained by some means, and those that do factorize are only partially known. 
Phenomenological analyses tend to either ignore power corrections or employ ad 
hoc estimates. On the other hand, a recent analysis (46) of power corrections 
within the framework of light-cone QCD sum rules (LCSR) found the possibility 
of large effects, of 0(10%) at the amplitude level. Importantly, some of the power 
corrections are of order A 2 /(4m 2 ), rather than A/m&. All this is sufficient to cast 
into doubt the cleanness of the optimized observables. 

In this paper, we revisit the issue of long-distance corrections, employing a 
systematic description in terms of helicity (as opposed to transversity) ampli- 



tudes. The point is that, when the treatment of 46 is reformulated in terms of 
helicity amplitudes, one can show a systematic suppression of the contributions of 
four-quark operators involving charm quarks to one (out of three) long-distance- 
sensitive helicity amplitudes, to be termed Hy. This suppression in A/m^ holds 
even at the level of A 2 /(4m 2 ) power corrections. A suppression is likewise ob- 
served for the contributions of the chromo magnetic penguin operator Q$ g to Hy. 
Moreover, if one adopts a resonance dominance model of long-distance contribu- 
tions at low q 2 , one finds again a suppression of the amplitude in question. This 
eliminates the leading candidate of large "duality-violating" corrections to the 
QCDF results (and power-correction estimates) for the latter contributions. 

The phenomenological relevance is as follows: Of the so-called "clean" observ- 
ables, there are precisely two which vanish in the limit Hy — > 0. These emerge as 
truly theoretically clean null tests of the Standard model, even with a conserva- 
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tive treatment of contributions from the hadronic weak Hamiltonian. Our results 
also show that this remains so at the q 2 — > end of the dilepton invariant mass 
spectrum, which translates into good sensitivity to the "right-handed" magnetic 



Wilson coefficient C 7 , as has been suggested earlier 72 , even when fully taking 



account of hadronic effects. For the other observables, however, the new physics 
sensitivity becomes more limited. 

The remainder of this paper is organised as follows. In Section [2] we set up 
our notation and review (and simplify) the decomposition of the decay ampli- 
tude and angular distribution, and "clean" observables, in helicity amplitudes. 
Section [3] contains a detailed discussion of sources of hadronic uncertainties both 
in the factorizable and non-factorizable contributions to the helicity amplitudes, 
and establishes the suppression of Hy. We also employ a parameterization of 
form factors at low q 2 which transparently separates the constraints from kine- 
matics and the heavy-quark limit from the issue of modelling power corrections. 
Section [4] comprises a detailed phenomenology of the "clean" observables, with 
particular attention to the low end of the low-g 2 region, which has traditionally 
been cut off ad hoc at q 2 = 1 GeV 2 . We find that the observables Pi and P^ p (in 
the notation of [57]) in B — > K*fi + fi~ stand out as theoretically cleanest, trans- 
lating to very good sensivity to right-handed currents, via the Wilson coefficients 
C 7 , C' g , and C' lQ . Specifically, we assess the (theoretical) sensitivity to the real and 
imaginary parts of C 7 to be on the order of 10% and 1%, respectively, with the 
sensitivity coming entirely from the region q 2 < 3 GeV 2 , and dominated by the 
g 2 -interval [0.1,2] GeV 2 . We also comment on the electronic mode, which shows 
a theoretical sensitivity to C' 7 very similar to the muonic mode. Throughout we 
take into account both the small but nonzero values of right-handed Wilson co- 
efficients in the SM and the effect of a nonzero muon mass, and show that two 
known algebraic relations in the massless case can be modified such that they 
hold, to excellent accuracy, in the presence of a finite muon mass all the way 
down to the kinematic end point. Section [5] contains our conclusions. 



2 Amplitudes and kinematic distribution 
2.1 Weak Hamiltonian 

The process B(p) — > M(k)£ + £~ , where M is a charmless final state (not necessar- 
ily a single meson), is mediated by the AB = 1 weak effective Hamiltonian, which 
is a sum of hadronic and semileptonic parts (where "semileptonic" is understood 
to include the magnetic penguin terms), 

T-icS = Tics' 1 + ^eflh (4) 
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with 



l/had 
it, 



s/2 



p=u,c 



8fl 



j=3.. .6 



cff 



4G f 

"71 



A, 



C7Q77 + C 7 Q' 7 + CgQgy + CgQgv + CioQioa + ClO^lOA 



+ C' S Q' S + CpQp + C'pQ'p + CtQt + C' T Q' T 



(5) 



(6) 



The operators p are given in 73 , the are defined as 

e 



m b sa, lu P R F^b. 



16vr 2 

Q 9V = ^(s 7 P i& )(^) 



Qt 



4vr v 
4rr m w 



aP R b)(!l), 
sa^P R b){h^P R l) 



QlOA 

Qp 



9s 



167T 2 

^em 



m h sa, v P R G^b. 



4tt ^ /M 



(^Pl6)(/7V0a, 



(sP R 6)(l7 5 0, 
4tt mvi/ 



(7) 

and the primed operators Q\ are obtained from these by P R Pl, Pl — Pr in 
the quark bilinears. g s (e) denotes the strong (electromagnetic) coupling constant 
coming from the covariant derivative = d^ l + ieQfA fl + ig s T A A A {Qf = —1 for 
the leptons), a em = e 2 /{Ait) and m& the 6-quark mass defined in the MS scheme. 

The contribution of the semileptonic Hamiltonian Hf s to the decay amplitude 
factorizes (in the "naive" sense) into a sum of products of hadronic and leptonic 
currents, 



where 



l^cfrl^O — Ly avu+L^aAn+Ls as+Lp ap+Lj, L aTL,n+Lj, R aTR. 



(e+e-\Iri\o), 

L s = {£+r\U\0), 
i 



re 

TL 



{t + r\qJa^P L m, 



L» A = (i + r\h^ 5 i\o) 

L P =(£ + r\h 5 l\0), 
i 



(8) 



(9) 



n TR 



q v la»"P R l\0), 



and we have made use of the relation 

(scr, u P R{L) b)(la^P R{L) s) = \{sq v a^P R{L) b){lq p a^P R{L) l), 



(10) 



where q = p — k is the dilepton four-momentumj^] The hadronic currents ay^, . . . 
are expressed in terms of form factors and Wilson coefficients, and enter the 
helicity amplitudes given below. 



1 Equation (10) holds for arbitrary time-like four- vector q^. 
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The hadronic Hamiltonian 7i e s requires in addition two insertions of the elec- 
tromagnetic current (one hadronic and one leptonic) to mediate the semileptonic 

decay, 

^(had) = _f_ Jd*xe- ig - x {£n-\f™> lept (x)\0) Jd A ye iq - v {M\T{f m ^»{y)U^{Q)}\B) 
e 2 

HI) 

where j' em > had >A« = ^ e q q^q. Hence, while this contribution does not naively 



factorize, it can be absorbed into ay^ in (|8j). Before discussing the amplitudes 
in more detail, we comment on the approximations implicit in and some conse- 



quences of rt8ft, (11 ) 



• The semileptonic weak Hamiltonian is the most general one up to dimen- 
sion six and can accomodate arbitrary new physics with a heavy mass scale. 
This includes all the standard scenarios, such as supersymmetry, extra di- 
mensions and little Higgs. In the Standard Model, CV, Cg and Cio are 
sizable, C 7 is suppressed by m s /mb, and the remaining Wilson coefficients 
are negligible. 

• The hadronic weak Hamiltonian is the Standard Model one, neglecting the 
small electroweak penguin terms. Beyond the Standard Model, there is 
a large number of extra operators; however unless new physics effects are 
dramatic their impact (through a^ ad ) will be very small and we will ignore 
them below. Such scenarios are also constrained by hadronic B decay data. 

• We work to leading order in the electromagnetic coupling, but all formulae 
so far are exact in the strong coupling, with non-factorizable effects confined 
to a^ d . 

• The leptonic currents can be decomposed into spin-0 and spin-1 terms (L v , 
L^) or are pure spin-1 objects (L^ L , Lt^ R ). It follows that the dilepton can 
only be created in a spin-0 or spin-1 state. Angular momentum conservation 
then implies that A is also the helicity of M, which is thus constrained to 
the values ±1 or even if M has spin greater than onej^] 

2 This statement is exact, rather than a consequence of naive factorization, following from 
the well-known fact that a particle's orbital angular momentum does not contribute to its 
helicity. If M is a multiparticle state, eg Kir, we mean by "spin" the total angular momentum 
of M in its cm frame and by "helicity" the projection of the M angular momentum onto the 
total M momentum in the B rest frame. 
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2.2 Helicity amplitudes and helicity form factors 

We now carry out the decomposition of the leptonic currents in spins and helic- 
ities. The resulting coefficients give the well-known "helicity amplitudes" 



16 



A.l 



This is easily achieved 39 through the completeness relation (see Appendix 
for our conventions) 

= ^(MK(1,A). (12) 

A=±1,0 

Here e(l, A), A = ±1, 0, denotes a (spin-1) helicity triplet of polarisation 4-vectors 
for a vector particle of four-momentum g M and mass ■s/cp, and e£* = g M / yq^. We 
may picture the latter as the "time-like" polarization four-vector of an auxiliary 
virtual gauge boson of mass \/q^, but the decomposition works independently 
of the origin of the weak Hamiltonian, and also for the tensorial currents. The 
result is 



A = - £v(X)H v (\)- £a(A)#a(A) +L s H s +L P H P 

A=±1,0 A=±1,Q 

- £tl(X)H T l(\) - J2 C TR {X)H TR (X), 



(13) 



A=±1,0 A=±1,0 



where 



C V (X) 




H V (X) 


= <(A)<, 


£a(X) 


= eM&k, 


H A {\) 


= e*(A)<, 


Ctl(X) 


= e fJL (\)L l TL , 


Htl(X) 


= <(A)< L , 


Ctr(X) 


= e ^)L TR , 




= e*(A)a Ti? , 




= L S , 


H s 


= as 


c P 






2m e a 


= L P , 


Hp 


= a P + q^a'l 
q 2 



(14) 



We have made use of the fact that all leptonic currents except for are con- 
served, so e t ^ contracts to zero with them. Moreover, the axial current obeys 
q^L^ = 2miLp, which allowed us to absorb the spin-zero axial vector amplitude 
into Hp (39]g 

The helicity amplitudes Hy, Ha, Hp, Hs are related to the "standard" helicity 
amplitudes fl8l|39] as follows, 



H XL/R = t V r f^(H v (\) T H A (\)), A t = i^-y/fH P , A s = -i y/fHs, 

(15) 



3 We do not distinguish between the lepton mass mi and the lepton field mass parameter, 
as we will work to leading order in the electromagnetic coupling. 
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where / is a normalization factor, which for M = K* and the conventions of 39 



is equal to F defined in Section 2.3 below. The helicity amplitudes H±\^r) are 
often expressed in terms of transversity amplitudes, 



.4 



\\L{R) 



V2 



±L(R) 



V2 



+i,L(R) 



H-iMR))- (16) 



However, we will work with helicity amplitudes throughout this paper, for reasons 
to become clear below. Explicitly, we have 



H v {\) 

H A (\) 
Htr(^) 

Htl(X) 

H s 
H P 



-iN\C 9 V LX + C' 9 V RX + 

iN(C 10 V L x + C' 10 V RX ), 
ArhbTriB 



m 



B 



2 m b 



{CjTl\ + C 7 Tp\) 



-iN 



CtT, 



m w v 9 



. Am b m B ~ 
-iN ^C T T RX , 



iN 
iN 



m Ws /q 
rhb 



m w 
rh b 



(CsSl + C' s Sr), 
{CrSl + C'pSp) 



m w 
2 m(fh b 



Cio ( Sl -Sr 



C 10 ( Sp 



m s 5 
m b 



(17) 
(18) 

(19) 
(20) 
(21) 

(22) 



where 



is a normalisation factor, 



N 



In 



AGfttib e 2 
y/2 16tt : 



rn 



(23) 



B 



contains the contribution from the hadronic hamiltonian, i.e. all non-factorizable 
effects, and we have defined helicity form factors 



- im B V L(R) x(q 2 ) 
m 2 B f L{R) x(q 2 ) 
im B S L ( R )(q 2 ) 



(M(\)\sf(\)P L{R) b\B), 
e*»(\)q»(M(\)\sa^P R{L) b\B} 7 
(M(X = 0)\sP R(L) b\B). 



(24) 
(25) 
(26) 



These expressions are still general enough to describe an arbitrary charmless final 
state M. Concretely, for a two-spinless-meson final state, not necessarily origi- 
nating from a resonance, the form factors will carry dependence on the dimeson 
invariant mass k 2 and its angular momentum L, in addition to the dilepton in- 
variant mass q 2 . 
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Note that parity invariance of strong interactions implies the relations 



V LX 
Tlx 
Sl 



-v( 

-T)( 
- V (-l 



\rv R ^x = v x , 

D^T R x = T\, 



Sn = S, 



(27) 
(28) 
(29) 



where s and rj are (respectively) the angular momentum and intrinsic parity of 
M. For a resonance, its spin s replaces L. Hence there are seven independent 
helicity form factors for spin > 1 and three for spin (when A = 0). Helicity 
form factors have previously been used in the literature as a technical vehicle 
in constraining form factors from unitarity 74 . As we will explain in detail 



below, helicity form factors are also preferable over the standard basis for form 
factors in weak decays: Not only do they simplify the expressions, but some of 
them are systematically suppressed, which can and should be exploited to reduce 
important sources of uncertainty. 

We also find it convenient to define rescaled helicity-0 form factors as 



Vo(q 2 ) 
To(q 2 ) 

S(q 2 ) 



2m B *Jq 2 ~ r t 2^ 



AV2 
2m B 



-w 

:% 2 ), 



/g2 A l/2 

2m B (m b + m s 
W 2 



~S(q 2 



(30) 



where A = 4m 2 3 |/;;| 2 (k is the 3- momentum of the recoiling meson in the B rest 
frame), and also define V±i(q 2 ) = V±i(q 2 ), T±i(q 2 ) = T±i(q 2 ). The helicity form 
factors can be expressed in terms of the traditional form factors. For a vector, 
we then have (conventions for polarisation vectors and form factors in Appendix 
Al) 



V±(q 2 ) 
V (q 2 ) 

T±(q 2 ) 

Uq 2 ) 
S(q 2 ) 



m B J m B (m B + m v 

1 



V{q 2 



2my\ 1 / 2 {m B + m v ) 



[{m B + m y ) 2 (m| - q 2 - m^A^q 2 ) - XA 2 {q 2 )} 



'"1: ~ m V T( 2, 2 
T2{<1 ] 2mJ Tliq } ' 



t 2 rrp 2 

2mT 



1 B 



1 13 



m B 



2m v X 1 / 2 
Ao(q 2 ), 



[m 



B 



3m 2 v - q 2 )T 2 {q 2 ) - 



A 



We also have Vrx = — V- X , T t 



R\ 



(m 2 B - m 2 v ) 
-T-x, S B = —Sl- 



Uq 2 ) 



(31) 
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For a pseudoscalar, we have 



V (q 2 ) = zf + (q 2 ), (32) 

To(d 2 ) = V ^ ^ (g 2 ), ^ 
{m B + mp) 

1 _|_ m s 2 2 

S(f) = l JL mB ~ 1 ™ M htf)- (34) 

m b 



In this case, V R0 = V Q , T R0 = T , S R = S. 



2.3 Kinematic distribution 

We now consider the process 

B{p) -> V{k)[-+ KfaMkt)] r( qi )£ + (q 2 ), (35) 

i.e. decays to a vector decaying further into two pseudoscalars (for definiteness, 
a kaon and a pion, with K = K° or K~ , tt = tt + or 7T°), where all four final-state 
particles carry definite four-momenta. These are the states in which detection is 
made. Following [32], we define angles 9k, 9i, 4> a s follows. We first define, in the 
B rest frame, 

Pi X Pl+ PK X Pn . PR + Ptt 

ei = | r, e A: = i r> 2 = i , r- ( 36 ) 

Then define 0, in the interval [0,27r], through sin0 = (ej x e^) • 2 and cos0 = 
eyex- Moreover, 8i is defined as the angle between the direction of flight of the B 
and the l~ in the dilepton rest frame and 9k as the angle between the direction of 
motion of the B and the K in the dimeson (K*) rest frame, both in the interval 
[0, 7r). For a B decay, we define the angles in the same way, in particular Q\ is 
the angle between the t~ (rather than the £ + ) and the B. This convention agrees 
with 1 39 1 and leads to simple expressions for untagged observables. 

Next, we assume a resonant decay through an on-shell vector meson. (This 
means we are making a narrow-width approximation^ We should then make 
the replacement 

\K*\ A) — ► VbJ <m K Y?{9,<l> K )\9 K]( t) K ), (37) 

where 9k is the angle between the +z direction and the K direction in the K* 
cm frame and <$k is the angle between the x axis and the projection of the 



4 Off-resonance effects vanish in the limit of an infinitely narrow K* . They can be included 
in the framework by introducing dependence on the hadronic final state invariant mass k 2 and 



55 



. In particular, off-resonant L = 



total angular momentum L in the helicity amplitudes 
(S'-wave) contributions have been recently studied in 
angular coefficients, but do not impact on those that involve only A = ±1 amplitudes. 



5N 



61 62 . They modify some of the 
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former onto the xy plane, and b = BF(K* — > Kir) m 1. (Except for the zero 
point of the angle (fix, this is entirely fixed by conservation of probability and 
of angular momentum and is independent of the details of the K* decay vertex. 
See also |55|.) Squaring the amplitude and summing over lepton spins the fully 
differential decay rate is obtained as 



9 



(38) 



dq 2 d(cos6i)d(cos6 k )d(f) 32 n 

x (i* sin 2 9 k + I\ cos 2 9 k + (JJ sin 2 9 k + l c 2 cos 2 6 k ) cos 29 1 

+I3 sin 2 6k sin 2 61 cos 20 + 7 4 sin 29 k sin 29\ cos 
+I5 sin 29k sin 61 cos + (Jg sin 2 6 k + Jg cos 2 cos 61 

+I7 sin 2$^ sin sin + lg sin 2^ sin 26 1 / sin + J 9 sin 2 0& sin 2 sin I 

The angular coefficients /j are functions exclusively of q 2 . They can be expressed 
in terms of the helicity or transversity amplitudes Eqs. (17)-(22) as 

n = ^{^(i<i 2 + i^i 2 ) + i^pI 2 + ^(i<i 2 -m 2 )+/3 2 ™ 2 }, 

Hy-\ 2 + \Hy\ 2 + (V^A)) + T ^ (|tf+| 2 + \H V \ 2 -(V^ A)) } 



J? 



/3 2 +2 



1 2 



1 2 



-F^(\H^+\H^), 
F^{m\ 2 +\Hv\*) + (V^A), 



h = --Re[H+(Hyy]+(V^A), 

h = F^Re[(H-+H+)(H v y]+(V^Al 

h = FURe[(H--H+)(H° A )*]+(V^A) 

I* = FPRe[H v (H A y-H+(H+y], 
It 



(3m e 



Re[H* s (H+ + H- 



2F^Re [H* S H V ] 



h = F^lra[(Hi + H A )(H v y + {V A)] 

Is = F^Im[(H v -Hy-)(H v y]+(V^A), 
h = F^lm[H+(H v y]+(V^A), 



Im [H*(Hy - HJ 



(39) 
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where 



AV W BF(K* -> A'.), 4 '"' 



3x2 5 vr 3 m 3 R v n \l q 2 



We have omitted the terms involving the tensor amplitudes Htl an d Htr, which 
will be considered elsewhere. The phenomenology in Section [4] will be incomplete 
in this sense, but as explained there, we expect the effect on the considered 
observables to be very small. 

The analysis of the CP-partner decay B — > K* // + /i~ gives rise to an equiva- 



lent distribution, d^'T/(dq 2 d(cos9i)d(cos 9k)d(fi), which is obtained from (38) by 
replacing 

^l«(c),2s(c),3,4,7 ~~ * hs(c), 2s(c), 3,4,7 j ^5,6,8,9 — ► ~h,6,8,9: (40) 

when one uses the angles defined as in the B decays with K~ — > K + fl8]. In 
the later equation the Jj's are equal to the Ij's but with all the weak phases 
conjugated. 



2.4 "Clean" observables, helicity hierarchies, and new physics 



The 12 angular coefficients, together with their CP-conjugates, provide complete 
information about the decay distribution. In practice, however, certain combi- 
nations of these observables a re more u seful. Two important examples are the 
forward-backward asymmetry 
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15 



17 



(2) 



and the transverse asymmetry A T 



32 



A FB {q 2 ) = 



,(2) = 



£ ~ J-il d(cos O^T' /(dq 2 d cos 6i 



/_! d{cos 9i)d 2 V j [dq 2 d cos 9\ 



3E| 
4T' 



2S 



(41) 
(42) 



2.s 



where Ej ee (/j + Jj)/2 and c?r' = (dr + dT)/2 denote CP-averages. These 
have two potential advantages. First, taking ratios leads to cancellations of form 
factor uncertainties between the numerator and denominator. In particular, Pi 
only involves the amplitudes Hy and H^. In the factorizable approximation and 
neglecting the small Wilson coefficients C 7 , C' g , and C[ Q , this leaves dependence 
only on T_, V-, T+, and V + . For energetic Ey ~ (small q 2 ), one has the 
heavy-quark/large-energy relations T_ = V-,T + = V + = (see below), which 
are broken only by (calculable) a s -corrections and (incalculable) A/m^ power 
corrections. Neglecting these, too, the single remaining form factor then cancels 
out between the numerator and denominator. In the case of Afb, the cancellation 
is incomplete, as the denominator involves further amplitudes. However, the 
location of the zero-crossing q 2 is determined entirely by the numerator, and 
(under the same approximations) is free of form factor uncertainties. For these 
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reasons, Pi, q%, and similar observables are often termed "clean". An optimized 



set has been recently defined in 57 and will be studied in the phenomenological 
part below. 

A second point is that Pi actually vanishes under the stated approximations, 
as a consequence of all terms being proportional to either V + or T + . Hence, it is 
an approximate null-test of the Standard Model, and a probe of any new physics 
that generates the Wilson coefficients C' 7 , C' g , or C[ Q . The same is true of J 9 and 
certain combinations constructed from it. 

Clearly, the actual theoretical cleanness of the observables will depend on 
the size of the radiative and power corrections and non-factorizable effects. The 
following section is devoted to a thorough study of these effects, and their impact 
on the "wrong-helicity" amplitudes Hy and H\ in particular. We will show that, 
under very conservative assumptions, Hy and H\ remain suppressed, such that 
the clean character of I3 and 1$ as null tests, but not of other observables, is 
preserved by non-factorizable and power corrections. 

Finally, let us recall that the radiative decay B — > V'j is described in terms 
of a subset of the amplitudes for B — > V£ + £~ . The precise relation is (A = ±1) 

a 2 

A(B V(X)j(X)) = lim -H v (q 2 = 0; A) 

g 2 -»0 e 



iNm 2 B 



^(C 7 T A (0) - C;T_ A )(0) - 16n 2 h x (q 2 = 0) 
m B 

(43) 



3 Helicity amplitudes: anatomy, hierarchies, and 
hadronic uncertainties 

The helicity amplitudes governing the observables involve form factors and the 
nonlocal objects h\, all of which carry hadronic uncertainties, limiting the sensi- 
tivity of rare B decays to new physics. However, hadronic uncertainties can be 
constrained by means of the equations of motion, the V — A structure of the weak 
hamiltonian, and an expansion in A/mb (QCD factorization). Our main point is 
that this results in the suppression of entire helicity amplitudes, including non- 
factorizable effects, such that the discussion is indeed best framed in terms of 
helicity (rather than transversity) amplitudes and helicity form factors. We first 
translate what is known about the form factors to the helicity basis, including 



the fact that the heavy-quark limit implies the suppression of two of them 20 
We next survey how this bears out in various theoretical approaches to form fac- 
tor determinations, concluding with a brief argument for the suppression of the 
positive-helicity form factors in the framework of light-cone sum rules, at the level 
of the correlation function. We then show that the V — A structure also implies 
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suppression of the "charm-loop" contribution to the nonlocal positive-helicity am- 



plitude h +1 , building on a method introduced in 46 . In addition, we show that 



the same conclusion applies to hadronic resonance models for the "light-quark" 
contributions to h\, once known experimental facts about the helicity structure 
of B — >■ VV are incorporated (which can be theoretically understood on the same 
basis). 



3.1 Form factors 



The B — > M form factors are nonperturbative objects. In the following, we re- 
strict ourselves to the B — > V case. First-principles lattice-QCD computations 



are becoming available 75 , 76 , although they will be restricted for the foreseeable 
future to the region of slow-moving V (high q 2 ). A state-of-the-art method of 
obtaining form factors at low q 2 is given by QCD sum rules on the light cone 
(see 69 77]). This involves, unfortunately, certain irreducible systematic un- 
certainties which are difficult to quantify. Sum rules are also useful in guiding 



extrapolations of high-g lattice-QCD results 74 



3.1.1 Theoretical constraints on form factors at low q 2 

The form factors fulfil two exact relations that in the helicity basis take the form 



S(q 2 



0) 
0) 



0, 

V o (0). 



(44) 
(45) 
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71 



At large recoil, i.e. small q 2 , one has further relations which hold up to correc- 
tions of 0(A/rrib) but to all orders in a s . As a result, the seven form factors are 
given, at leading power in A/m^ and A/E (where E = Ey is itself of order 
for low q 2 ), in terms of only two independent soft form factors 
with radiative corrections systematically calculable in QCDF 
bative expansion in a s . These corrections also involve nonperturbative objects 
such as decay constants and light-cone distribution amplitudes (LCDAs) of the 
initial and final mesons. The factorization properties and calculation of radiative 
corrections become particularly transparent when formulated as a matching of 



£l and 
as a pertur- 



QCD to soft-collinear effective theory (SCET) [78-81 . Corrections at 0(A/rrib) 
violate factorization and need to be modeled, or estimated or constrained by 
another method. 

The soft form factors can be chosen to coincide with two "physical", i.e. 
full-QCD form factors, which makes them well-defined to all orders in A/ttv 
Appropriate choices are V and A , for £j_ and £y respectively 71 , given that they 
are matrix elements of conserved currents and, as such, free from renormalization 
scale ambiguities. For we find it convenient to use the transversity form factor 
Ti instead. In this case, however, the transversal soft form factor depends on the 



14 



factorization scale /i as (to LL accuracy) 



Ti(g 2 ,/x) =£ ± (g 2 ,/i) = £±(q 2 ,m b ) 



Q s (/i) 
a s (m b ) 



4/23 



(46) 



Setting 



S(q 2 ) 



t±(q 2 ,v), 
— 5, fe 2 ) 



with £ ~ (mfj — q 2 )/ (2ms), the symmetry relations in Ref. 71 can be expressed 
in the helicity basis as 



6 





2E 








m B 


T+ 


= o, 




E 


r«» 




my 




2E 








m B 


v + 


= o, 




E 


V 




m y 



h 1 + 



a i + 



47T 



rnz 



In -4 - 2 + 4L 



en i + 



4tt 

47T 



In 



+ L 



mi 



47T 



ATn 



47T 



Ay, 



[-2 + 2L] 



47T 



AVq, 



(47) 



with = 4/3 and L = —2E/(mB — 2E) \n(2E /ms)- These expressions hold up 
to higher-order corrections in a s , which augment the a s terms shown, and power 
corrections. 

The a s -contributions multiplying the soft form factors come from the hard- 
vertex corrections, and the remaining ones originate from hard scattering with 



the spectator quark 71 



AT 

AV 
AV 



m 



B 



AE 2 

-,2 



AFii 



q 



AE 2 
2 - 



AFn 



(48) 
(49) 
(50) 



where AFj_ and AFy involve (finite) convolutions of hard-scattering kernels with 
light-cone distribution amplitudes and can be found in Ref. 71 . Higher-order 
corrections in a s do not change the structure or infrared safety of (47), i.e. fac- 
torization can be proven to all orders 82 . In particular, the vanishing of T + and 
V + at leading power is an all-orders result 20 , 83 , which looks this simple only 
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in the helicity basis. The 0(a 2 ) contributions have been calculated in 83-86 
and their numerical impact was found to be small, mainly reducing the residual 
(unphysical) scale dependence. 

Thus, in the heavy-quark/large-recoil limit, the form factors T + and V+ ex 
actly vanish. Combining (44) and (47), we have at low q 2 



T + {q 2 ) = 0(q 2 /m 2 B )xO(A/m b ), (51) 
V + (q 2 ) = 0(A/m b ). (52) 

On the other hand, T and V are not suppressed, and are independent of any 
hadronic information related to the transversal polarizations of the vector meson. 
Notice also that, with our choice of soft form factors, the vector form factor V, and 
hence VI, has a purely residual (higher-order) scale dependence fi at any given 
order of perturbation theory, from the factorization into the scale-dependent £j_ 
and a scale-dependent perturbative factor. One can explicitly check that this 
produces a relative change in the form factor of no more than a 1.5% in the range 
ra&/2 < fi < 2m b . 

3.1.2 Numerical values of the B — > K* soft form factors 

Although the symmetry relations reduce (at leading power) the number of in- 
dependent non-perturbative functions, for a quantitative treatment one still has 
to compute the soft form factors by a nonperturbative method and estimate (or 
calculate) the power corrections. Sum rules formulated on the light cone (LCSR) 
are customarily used in exclusive B decays to obtain numerical values of the 
form factors in the large- recoil domain |46j[69j[77]. Other approaches that have 
been used to calculate the form factors in this regime include local QCD sum rules 
(QCDSR) [87] and (truncated) Dyson-Schwinger equations (DSE) [88). We list in 
Table [T] the results on the B — > K* form factors at q 2 = and in the transversity 
basis for the different calculations considered in this paper. The central values 
of most of the form factors are quite similar, except some (prominently Aq) for 
which different methods disagree. 

A well known issue between the generic results of LCSR calculations (e.g. 
those of Ref. [77]), QCD factorization, and the SM value of the Wilson coefficient 
CV is that they lead to a branching fraction of the decay B — > K*^ that is larger 
than the experimental value. Given that CV is constrained to be close to its SM 
value by the inclusive B — > X s ^ decay rate (at least when assuming C' 7 = and 
Cj real), it is often assumed that this discrepancy is due to a systematic error in 
the LCSR model which produces an overestimation in the value of the relevant 



form factor Ti at q = 23 . A possible solution to this problem is to re-scale 



the form factors such that Ti(0), in combination with the SM value for CV, lead 



to the experimental branching fraction of the radiative decay 39 . Although this 



procedure discards in part the sensitivity of the B — > K*£ + £ decay rate to new 
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Table 1: Values of the B — > K* QCD form factors at q 2 = in the conventional 
basis and in the different approaches considered in our study. The uncertainties 
are those quoted or suggested in the corresponding references (a 15% relative 
error for the DSE results). 





LCSR 


77 


lcsr pe) 


QCDSR [87] 


DSE M 


V{0) 

4)(o) 
Mo) 

A 2 (0) 
Ti(0) 
Ts(0) 


0.411(33) 
0.374(34) 
0.292(28) 
0.259(27) 
0.333(28) 
0.202(18) 


0.36+^ 

0-25±g:ig 

23+ ' 19 
u - zo -o.io 

S1 +ai8 

u - oi -0.10 

22 +ai7 


0.47(3) 
0.30(3) 
0.37(3) 
0.40(3) 
0.38(6) 
0.3 


0.37(6) 
0.25(4) 
0.28(4) 
0.30(5) 
0.30(5) 
0.27(4) 



physics, it does not affect any physical information extracted from asymmetries 
or ratios, which are, indeed, the most interesting observables in the semileptonic 
decay [39~] . 

Using the SM value for C 7 and B(B ->■ K*j) = (4.33 ± 0.15) x lO" 5 , one 
obtains 



C±(0) =7i(0) = 0.277(13), 



(53) 



where the error is obtained adding linearly the experimental and the theoretical 
uncertainties, which come from the SM parameters and the hadronic parame- 
ters entering the non-factorizable contributions to the amplitude (shown in Ap- 
pendix A.2). Re-scaling the values of the form factors in Table [l] and taking the 
average of the results given by the different models, we obtain 



£n(o) 



2m 



K* 



m B 



S(0) = 0.09(2), 



(54) 



where we have estimated the uncertainty to be such that it includes the two 
extreme values coming from the LCSR 77 and DSE 88 . For the q 2 dependence 



we choose the pole forms introduced in 
arguments, 



23 based on pure heavy-quark-limit 



Z±(q 2 )=U0) 



1 — q 2 /m 2 3 



£n(9 2 ) = £n(o) 



1 — q 2 /m 2 3 



(55) 



In Fig. [T] we plot the form factors T\ and S in the large-recoil region used in 
this work as a reference for the symmetry relations in Eqs. (47) with the error 



bands produced by the uncertainties on £j_(0) and £||(0). These are compared with 



re-scaled central values of the results reported using LCSR 46[ 77 
and DSE 



QCDSR 87 
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Figure 1: Form factors T\ and S in the large-recoil region used in this work 
as a reference for the symmetry relations (thick red line) with the error bands 
produced by the uncertainties on £j_(0) and £||(0). These are compared with the 
re-scaled (see in the main text) central values of the results obtained in the LCSRs 
(blue dot-dashed [77] and black dotted [46]), in QCDSRs [87] (green short-dashed) 
and DSEs [88] (purple long-dashed). 



3.1.3 Power corrections to the large-recoil form factor relations 

The main factor limiting the utility of QCD factorization for the form factors 
(and elsewhere) are unknown power corrections 0(A/mb). Although one might 
naively expect such power corrections to be ~ (5 — 10)%, the possibility of having 
larger corrections cannot be dismissed. Let us refer to effects breaking the (QCD- 



corrected) symmetry relations in Eqs. (47) as power corrections, even though 



these also include the numerically unimportant perturbative a™, n > 2 contribu- 
tions. They govern the so-called factorizable power corrections in B — > K*£ + £~ 
decay. As they do not cancel out in the "clean" observables defined in [57], es- 
timating them is important to assess hadronic uncertainties. The conventional 
procedure in phenomenological analyses is to use the results of some technique 
that automatically includes power corrections (most commonly LCSR). However, 
in doing so, the systematic errors coming from the assumptions and approxima- 
tions implied by the particular approach are not transparent. A related issue is 
that the g 2 -dependence of the form factors is often "hard-coded" and not treated 
as an uncertainty. 

Instead, in this paper we parameterize power corrections to the form factors 
in a model-independent fashion. The uncertain parameters can then be estimated 
by various methods. For a given (helicity) form factor F, we parameterize the 



corrections to (47) as 




F pc- = a F + bF ^ + [l^) ;A 2 /m 6 2 ], (56) 



where and bp are dimensionless numbers of order (9(A/m&). Importantly, 
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Figure 2: Form factors in the helicity basis and in the large- recoil region used 
in this work (red thick and solid line). The inner (red) error band is produced 
by the uncertainties on £j_(0) and £||(0), while the outer (blue) includes also the 
estimated uncertainties on the factorizable power corrections. These curves are 
compared with the re-scaled central values of the results obtained in the LCSR 
(blue dot-dashed 77 and black dotted (46] lines), in QCDSR (87) (green short- 



dashed line) and DSE [88] (purple long-dashed line). For the sake of completeness 
we also include the curves (orange solid and thin line) obtained in the strict 
heavy-quark/large recoil limit, Eqs. (47). 



Eqs. (44) and (45) imply that ap + = and ay = 0. 

Ideally, the coefficients ap and bp should be estimated by a direct calcula- 
tion of helicity form factors, for example through LCSR (see discussion below). 
For lack of present availability of such results, we compare available results (in 
the conventional basis) of the LCSR, QCDSR and DSE approaches to give an 
estimate of the uncertainty given by the break-down of the symmetry relations. 
We estimate ranges for the nonvanishing ap and bp by taking the average devi- 
ations between the results given by the symmetry relations and those obtained 
in the different models considered in this paper. The resulting ranges are given 
in Table [2j Note that the of are all at the few-percent level, well in line with 
"naive" expectations about the size of power corrections. The (less important) 
bp coefficients are also around the 10-percent level or below, with the exception 
of the helicity-zero form factors Vq and, to smaller extent, To. In Fig. [2] we plot 
our ranges for the form factors, including power corrections described by (56), 
with the ranges for ap and bp given in Table |2j We also show the central values 



of the results obtained from LCSR 46 ,77 , QCDSR 87 and in DSE 88 



Our method may over- or underestimate the power corrections. For instance, 
one would obtain larger uncertainties by combining the LCSR uncertainties on the 
standard form factors entering V + and T + in quadrature (or linearly). Hence, our 
numbers should be considered preliminary until dedicated calculations appear. 



19 



Table 2: Bounds on the dimensionless constants dp and bp parameterizing the 
power corrections to the form factor heavy-quark relations. 







v + 


T_ 


T + 


v 


To 


I Q> | max 


0.027 


0.008 











0.050 


| b | max 


0.136 


0.042 


0.125 


0.043 


0.434 


0.206 



Nevertheless, it is reassuring that the nonperturbative approaches appear to be 
consistent with expectations from the heavy-quark limit. In fact, for the LCSR 
for the form factors the heavy-quark limit has been analyzed a long time ago 
in 89 , and the result is fully consistent with the structure later obtained in 71 
and expressed in the helicity basis through ({47]). 

To reduce the LCSR sum-rule uncertainty manifestly, one should derive sum 
rules directly for the helicity form factors]^] To demonstrate this, consider the 
correlation function 

G FX (q 2 ;p 2 )=i J dV^(^*W|T{e;( g ;A)(^6)[0]^(y)}|0), (57) 

where js = irribqjsb, q is a light quark field, and Tp a Dirac structure appearing 
in the definition of a given helicity form factor. Its hadronic representation, 

(K*(k,\)\e;(q;\)ST»b\B) f B ml 

Lrpx — 7, n r • • ■ i \°°) 

p z — m B nib 

obtained by inserting a complete set of hadronic states between the two cur- 
rents, contains the helicity form factor F\(q 2 ) oc (K*(k)\e*^(q; X)sT F b\B), with 
the ellipsis denoting contributions from higher resonances and continuum states 



with the quantum numbers of js- A sum rule is obtained from (58) by taking 
p 2 < m 2 B off-shell by an amount 0{1 GeV 2 ), evaluating Gp\ perturbatively in a 
light-cone expansion, and Borel-transforming both sides to suppress the hadronic 



states above the B. This is completely analogous to the standard sum rules 69 
and in fact simpler as the correlation function for a given V p and A directly gives 
the respective helicity form factor, with no need for Lorentz decomposition. The 
light-cone expansion of Gp\ results in convolutions of perturbative kernels with 
light-cone distribution amplitudes, which are organised in terms of increasing 
twist. Let us demonstrate the suppression of the positive-helicity form factors in 
this framework. At tree level, we have for T + 

»/dH/e-**T{^( g ;+)(gI^6)[0]ji(y)} = s(fci)e* (g; +)q^P R / ~[ 2 + ™» 2 g(fc 2 ; 
J \P ~ K i) — wib 

Now, to leading-twist accuracy we can replace s(k\) and q(k2) by collinear fields 



'Such direct calculations have been recently advocated, independently, by T. Feldmann 90 
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satisfying tfi_q = and sfi + = Oj^] As a consequence, we can insert a projector 
s — > s ^~f + . We then encounter in the Dirac structure a factor 

fi^ + ft*(q;+)iP R 

It is easy to verify that this vanishes as a consequence of 

This shows that the correlation function, and thus the form factor T + , is sup- 
pressed as 0(A/Ek*) = 0(A/mb) as a higher-twist effect. Independently, 

where at low g 2 <C m 2 ~ E 2 K * we have n + ■ q ~ m&, while n_ ■ g ~ g 2 /m 2 . Hence 
up to corrections of 0(g 2 ), T + vanishes as a consequence of 

f(q;+)P L ^ = 0. 



In combination, we obtain a behaviour that respects (51). In the analogous case 
of V + , one obtains the behaviour (52) in an analogous way, but no g 2 -suppression 



(for there is no $ factor in the definition of V + ). For a more quantitative treat- 
ment, one needs to go to higher- twist accuracy; this requires going beyond the 
collinear approximation for the external momenta and quark fields. In principle, 
the necessary ingredients are known up to twist-4. We leave such an analysis 
for future work. Note that the argument does not imply a suppression of the 
helicity-0 form factors. 

3.2 Non-Factorizable effects: helicity hierarchies for h\ 

We next turn to the contribution of the hadronic weak hamiltonian in rare 
semileptonic -B-decays. As explained above, these contributions enter only in 



ify(A), via h\ defined through (11) and (23). They do not naively f actor ize 
into form factors and leptonic currents, but involve the B to V matrix ele- 
ments of a T-product of the weak hamiltonian and an electromagnetic current 
(1/|T(j em ' had ^( 2 /)^ d (0))| J B). A systematic treatment based on QCD factoriza- 
tion exists [23|[29] , to leading power in an expansion in A/mj,. Contributions 
can be divided into "form factor correction", "annihilation", and "spectator- 
scattering" terms and are given in terms of (soft) form factors, hard-scattering 
kernels, and light-cone distribution amplitudes, as in the form factor relations 
discussed above. Power corrections do not, in general, factorize. Only a subset 



6 We are assuming the K* to move in the +z direction here, as throughout the present paper. 
Note that in the relevant literature on LCSR for form factors, the final-state hadron is often 
taken to be moving in — z direction instead. 
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of A/mb power corrections to B — > V£ + £ have been studied in the literature, in 



the context of isospin asymmetries 27,28 . Some of these contributions violate 



factorization and were modelled in a standard way 27,91 . 

Importantly, QCD factorization predicts that the suppression of the positive- 
helicity amplitudes continues to hold in the presence of non-factorizable terms 



23 , and it is also true for the known power corrections 27 28 . That is, h + still 
vanishes at this order, and this is again a direct consequence of the V —A structure 
of SM weak interactions. For h_, the known power corrections do not vanish, 
and result in an isospin asymmetry of about 10% at q 2 = 0, rapidly falling with 
increasing q 2 . For the isospin-averaged branching fraction, this reduces to about 
a 3% effect [from a replacement e u — td — > e u + . However, the modelling of 
(incalculable) power corrections is not necessarily accurate, and it is in particular 
important to investigate to what extent h + may be generated by them. 

In the following, we separate the hadronic weak hamiltonian into a charm 
part (with large CKM and Wilson coefficients), a gluonic part (Qsg, with a large 
CKM and Wilson coefficient), and a light-quark part (involving only operators 
which are suppressed by small CKM elements or by small Wilson coefficients), 
and investigate in each case the possible sizes of power corrections, in particular 
to h + , by methods complementary to QCDF. 

3.2.1 Charm loop 

Within the context of LCSR, a study of charm loop effects at low q 2 has been 



given recently by Khodjamirian et al 46 , and the analogous contributions to 



B — > K*^/ have been considered earlier in [34] . In |46], long-distance charm-loop 
effects are estimated to be sizable (and with large uncertainties); these effects 
correspond in part to power corrections in QCDF. Unfortunately, the results 
in 46 are only presented in a numerical form and only for transversity, not 
helicity, amplitudes, expressed through effective (amplitude-dependent) shifts of 
the Wilson coefficients Cg. Nevertheless, central values and uncertainties given 
there are suggestive of a strong correlation of the charm loop contributions to -An 



and A± and a suppression h + <^ h_. The computation in 34 provides directly 



results for h + and /i_ at q 2 = 0, both of them numerically smaller, implying also 



h_ much smaller than suggested by 46 



The aim of this section is to argue that a hierarchy h + -C h_, ho results, as 
far as the charm loop goes, from the light-cone dominance of the amplitude at 
q 2 <C m 2 B . To this end, let us recast the strategy of 46 in terms of helicity 
amplitudes, picking out the charm loop in h\, 

h x \ c5 = ~e;(X) I 'd*ye i ^(M|r{(cyc)(y)(Cjgf + CS02)(0)}|S>.(59) 
m B 6 J 



Next, 46 shows that the Fourier integral is dominated by the light-cone y ~ 0. 



A light-cone OPE is then performed. To leading order, this results in a local op- 
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erator whose matrix elements can be identified with the charm-loop contribution 
to the form factor term in QCDF (ie those charm-loop effects that do not involve 
the spectator quark). At the one-gluon level, one has the expression 

Mcc,ld = e^(X)(M(k,X)\d^\B), (60) 

where 

6, = J duI, paP {q,uj)7s L Y5(uj - % ^±l^G a ^b L , (61) 
with D the covariant derivative and InpaB given in [461. The nonlocal operator 



(61 ) is the first subleading term in an expansion in A 2 / (4m 2 — q 2 ), with terms in- 



volving two and more gluon fields contributing only at higher orders [46] . Eq. ( 60 ) 
hence provides an approximation to the long-distance charm-loop contributions. 
It can be further expanded in local operators, 



^ - ^w^)L/,r(^)>%. (62) 



in a. ■ D\ n ~ 



The result of [34] corresponds to keeping only the n = term, and evaluating its 
matrix element by means of a LCSR for a correlation function 

• [d*ye~^(K*\[d^(q)}(0)jUy)\0). (63) 



Ref. 34 argued the suppression of higher terms in the local OPE by a larger 
expansion parameter of order m j eA/(4m 2 ), which has been taken as (20 — 40)% 
and used to justify truncating the OPE after the leading term. This numerical 
value corresponds to taking A ~ 300—650 GeV (for MS quark masses), and should 
hold up to an 0(1) factor, which if large could in principle spoil the convergence 
of the OPE. More seriously, the power counting itself was obtained by appealing 
to inclusive B — > X s 7 decay, where similar matrix elements (B\b(q-D) n G a pYb\B) 
occur as part of power corrections to the charm loop [92j[93]. (r denotes a Dirac 
structure which is irrelevant to the present discussion.) There, the softness of 
the B meson constituents provides one power of A in the numerator, which can 



be seen via q ■ D D —iq ■ ka ~ m&A, where kc is the gluon momentum 92 



(The resulting 'suppression' factor is estimated as 0.6 in |93|.) However, with an 



energetic K* in the final state as in (63) the constituents have energies 0(mb), 



so n + ■ D D n + ■ k G ~ m b and a scaling m 2 /(4m 2 ) of the putative expansion 
parameter seems appropriate; at least, establishing a suppression requires a new 



argument. We therefore will not rely on the estimate of 34 in this paper. Ref. 46 
estimates instead the full nonlocal operator matrix element from a LCSR for a 
different correlation function 

(0\T{jr(y)O,(0)}\B), (64) 
where jff* = d^^s, which yields the matrix element in terms of 5-meson LCDAs. 
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Figure 3: Tree-level diagram for the light-cone OPE for the correlation function 



(65). The dashed line indicates a two-particle cut contributing to the (perturba- 



tive) spectral density. 



To show the suppression of h + , note that h± can be obtained directly from 
G hX (q 2 ; k 2 ) =-i f dV^(0|TK*(5; A)jf (y)e^(-z; X)O,(0)}\B). (65) 



To be precise, we take k = (k°, 0, 0, |k|), as well as q, in the {tz) plane. Note that 
for A = ± the polarisation 4-vectors are (with these conventions) independent of 
k, hence the rhs indeed defines a Lorentz- invariant function of k 2 and q 2 . (The 
formalism could, with appropriate care, be extended to A = 0.) The hadronic 
representation contains the desired matrix element, 

n l 2 ; 2\ fK*\\ TTT'K* 

G hX {q ; k ^ 



m 2 K , — k? 



(K*(k; X)\e^*(-z; A)O M (0)| J B) +continuum contributions. 

(66) 

Here k = (\/m 2 K * + k 2 , 0, 0, |k|) is the physical (on-shell) 4-momentum of the 

we take 



K* corresponding to the given q 2 . To obtain a LCSR, following 46 
k 2 ~ — lGeV 2 ~ —m^A (corresponding to a Borel parameter ~ ^m^K) and 
consider the light-cone OPE of Gh\. The leading (tree) diagram is shown in 
Fig. [3j Defining n + ■ q = itlb + I, and taking the 6-quark momentum to be 
Pb = rrib/2(n + + n_) + ri, we have / ~ A > and 

,2 



q 



2 m b + I 2 



k 

Q 



n 



_+ 
2 

n + 



+ 0(A), 
+ 0(A), 



(67) 
(68) 
(69) 



such that Q 2 ~ m&A, making the internal propagator, in QCDF/SCET terminol- 
ogy, a hard-collinear line. The operator product, to tree level, becomes 

-z;\)0,(q), (70) 
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where Pr has been inserted for convenience [note that 0^{q) contains a chiral 
projector]. Neglecting terms 0(A,m s ) in the propagator, this vanishes for A = +. 
Hence the nonlocal charm loop contribution to h + is (9(A/mb)-suppressed relative 
to h-, and altogether 

fc+i — = ° < 7i > 

relative to the leading-power amplitude Hy. Nevertheless, we will see below that 
it constitutes the dominant remaining source of theoretical uncertainty on several 
"clean" observables, thanks to the strong suppression of the helicity form factor 
T + . This shows at the same time the importance of taking it into account, and 
motivates further, more quantitative work on the nonlocal power-corrections to 
(mainly) the positive-helicity amplitudes. 

in- 



Unlike in the form factor case, the more traditional LCSR of Ref. 34 



volving the correlation function (63) with an on-shell K* rather than a B, does 
not lend itself to an analogous argument. Essentially, the reason is that that 
sum rule is given in terms of (most importantly) chiral-even twist-3 3-particle 
LCDAs, where, loosely speaking, the helicity of the meson is determined by that 
of the gluon field-strength tensor, whereas the helicities of the quark and anti- 
quark fields cancel out. Hence the presence of a chiral projector has no direct 
impact on the helicity of the K*. Of course the suppression should still be seen 
in the numerical result for the full (rather than lowest-order local OPE) result. 
However, higher-order terms in the local OPE are not known, nor do they seem 
to be suppressed, as argued above. 

3.2.2 Chromomagnetic penguin operator 

The operator Q$ g comes with a large CKM factor and its Wilson coefficient is 
not small. In QCDF, it enters exclusively through the hard spectator scattering 
contributions. These factorize at leading power, but at order A/m b , attempting 



to factorize them results in endpoint divergences 27 , implying a breakdown of 



factorization. Cutting off the end point divergence and modelling the end-point 



contribution as in 27 s 28 , h+ still vanishes at subleading power. However, on 
general grounds one would expect chirality-violating QCD (and strange mass) 
effects to generate a contribution to h + at order A/m b and m s /m b . It is possible 
to estimate the Q$ g contribution to h\ with a LCSR, similarly to the sum rule 
for the nonlocal operator M . For B — > K£ + £~ (where the dilepton helicity 



A = 0), such a calculation has been very recently performed in 94 , wherein a 
"soft" contribution is identified and expressed in terms of three-particle 5-meson 
light-cone distribution amplitudes. The operator M is replaced by the simpler 
(local) Qs g , but the required correlation function involves an extra insertion of 
the electromagnetic current anywhere on the quark line in Fig. [3j 

Here, we note that our above argument showing a helicity hierarchy of the 
long-distance charm loop contributions likewise can be applied to this sum rule 
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for the matrix element of Qg g . The operator Q 8g still provides a chiral projector, 
and the fermion line entering j K * in Fig. M is still "hard-collinear" , such that 



the first four factors in (70) are unchanged, even if the electromagnetic current 
insertion occurs on that line (if the strange quark mass is neglected). The result is 
again that the long-distance (soft) contribution /i+|q 8Si ld is suppressed by A/m^ 
or m s /mt> relative to /i-|q 8S ,ld- As the latter is already suppressed by a power 
of A/rrib relative to the leading-power amplitude Hy, the impact on Q Sg on h + 
should be negligible. (Note also that the effect of the "soft" Q 8g contributions in 
B — > K£ + £~ was found to be well below 1 % of the total hadronic contribution 



m 



94 . It is difficult to see how a much larger contribution could occur in the 



present case, even for the non-helicity-suppressed amplitude 7t_.) 
3.2.3 Light quarks and resonance structure 

The remaining contributions of the hadronic weak Hamiltonian to the decay 
amplitude coming from the QCD penguin operators and the double Cabibbo- 
suppressed current- current operators involving up quarks, 

a had,i q = J # xe -*r*(K*\T{fi*(x), n^ d ' l \0)}\B), (72) 
and are either doubly Cabibbo-suppressed or weighted by the small Wilson co- 



efficients C*3_6- Again, a systematic description exists within QCDF 23 , with a 
vanishing contribution to h + at leading power and a breakdown of factorization at 
subleading powers. Because of the multiple suppression factors, the contributions 
to Hy arising in this fashion are negligible. 

However, long-distance non-perturbative effects may manifest themselves partly 
as resonances or poles in the complex-g 2 plane, implying a resonance structure 
which we do not expect to be accounted for at any order in A/m&. Therefore 
we employ a hadronic description to estimate both power corrections and the 
possibility of large "duality-violating" effects in B — » K*£ + £~ observables. In 
order to do this, let us consider instead the object 

~had,i q = J d 4 xe- iq - x (K*\T{jp lci (x), nX(Q)}\B), (73) 

where we only keep the light-quark part of the electromagnetic current, relevant 
for resonance structure in the low-g 2 region (but revert to the full weak Hamil- 
tonian). Ideally, we would like to compute a^ ad ' lq taking into account the fact 
that pions and other light hadrons are the relevant degrees of freedom of QCD 
in this domain, in a systematic fashion as, for example, using chiral perturbation 



theory (xPT) [95,96 , together with any of the methods that extend its range 



of applicability up to the region of the light resonances 97 100 . In fact, this 



program is attainable for kaon decays in which the energies and masses are all 



small compared with the chiral symmetry breaking scale A x s# ~ 1 GeV 101 . In 
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f V B K' 




Figure 4: Graphical representation of the VMD model. The filled bulb represents 
the B — > V K* decay vertex, as obtained in QCD factorization, the solid bulb fy, 
as obtained from experiment, and the double lines resonance propagators, with 
the cross indicating the multi-particle dressing of the respective pole. 



case of the B- meson decays, one encounters heavier scales and, at present, it is 
not clear how to integrate them out in order to extract the long-distance effects 



in a model-independent way 102, 103 



In this paper we use a model to estimate the contribution of the light hadronic 
degrees of freedom in the low-g 2 region. We start by making a factorization 
approximation of the correlation function Eq. (73), using a basis of hadronic 
states |P(0)) and \P'(x)), 



~had,lq 



^ e -^x<°br^)i^^ ( 74 ) 



P,P' 



where the sums include further integrations for multi-particle states. We next 
assume that these sums are saturated by the lightest neutral vector resonances 
V = p(770), w(782) and 0(1020), i.e. vector meson dominance (VMD). This 
hypothesis has proven very fruitful in modelling the electromagnetic structure 
of light hadrons at low energies. It finds microscopic justification in the large 
iV c limit of QCD 104 and it has been successfully implemented to connect the 



short-range part of the low-energy interactions of pions with QCD 97,98]. (For 
a compilation of phenomenological applications of the model in the weak decays 
of mesons see Ref. 



105 



In the VMD, the first factor in the RHS of Eq. (|74j) is a 
semileptonic decay constant, fy, the second the vector-meson propagator and the 
third a B — > VK* decay amplitude. Finally, we (partially) take into account the 
effect of the continuum of multi-particle hadronic states by dressing the poles of 
the resonance by their (off-shell) width. All in all, the estimate for the hadronic 
contribution at low q 2 can be pictured as in Fig. |4j 

In order to carry out the computation, it is convenient to use an effective La- 
grangian containing fields which serve as interpolators for the vector resonances. 



We choose the anti-symmetric representation advocated in Refs. 97, 98 1 for appli- 



cations in xPT. Other Lagrangian formulations consistent with chiral symmetry 
and electromagnetic gauge invariance^Jare equivalent to this one, once consistency 

7 Notice that in a previous VMD analysis |l06| of the vector-meson contribution to the B — s- 
K*£ + £~ decay, electromagnetic gauge and non-gauge invariant Lagrangians were considered in 
the same footing and large differences between the two approaches have been reported at low 
q 2 . In this paper we work exclusively with approaches consistent with electromagnetic gauge 
symmetry (and QCD, as stated in the main text). 
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Figure 5: Relative contribution (see main text) of the light-quark component of 
the electromagnetic current to the helicity amplitudes HY and Hq. The dashed 
(blue) line is the result in QCD factorization while the solid (red) line is the result 
in the hadronic model used in this work. 



with QCD asymptotic behavior of 2-point spectral functions is demanded 98; . 
We address the reader to Appendix [B] for the details and conventions used in the 
model. 

As for the B — > VK* decay amplitude, it is natural, in the present context, 



to use the QCD factorization calculation reported in Ref. 107 . In fact, as al 



ready discussed in 23 , there is a one-to-one correspondence between a subclass of 
diagrams in the QCDF calculation of B — > K*£ + £~~ and of the diagrams appear- 
ing in the QCDF calculation for B — > VK*, where V is one of the light-quark 
resonances. For our model calculation, we employ the QCDF results for the 
B — > VK* amplitudes, while at the same time omitting the corresponding terms 
from the B — > K*£ + £~ amplitude. Next, the heavy-quark limit predicts a hier- 
archy of helicity amplitudes, Ho : H : H + = 1 : A/m^ : A 2 /ml, i.e. in relative 
terms H + /H_ ~ A /mi,. Although H_ and H + do not factorize in QCDF, as the 
hard-scattering contributions develop end-point singularities, and the hierarchy 
Ho ^> H^ is numerically removed by large penguin corrections (obviating the 
so-called polarization puzzle 108] ) , the hierarchy H + H^,Hq remains and is 



in good agreement with experiment in B — > K (f> decays, for which complete 



angular analyses exist [109 110 . An important consequence of this is that both 
in QCDF and in our model, the contribution to the B — > K*£ + £ "wrong-helicity" 
amplitude is suppressed as 



h h + a(i ' lq /h^' i(l = 0(A/m b ) 2 , ^ ad ' lq /^ ad ' lq = C(A/m fe ), (75) 

on top of the smallness already implied by the CKM factor and Wilson coeffi- 
cients. As a result, we expect the light-quark hadronic pollution of observables 
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sensitive to the chirality-flipped operator Q' 7 to be negligible, and preclude in par- 
ticular large "duality-violating" corrections to the leading-power QCDF results 
for B K*£ + £~ . Thus, in the phenonomenological section we set /i^ ad ' lq = 0. 

In Fig. [5] we compare the model calculation to that obtained from the corre- 
sponding subclass of QCDF diagrams. More precisely, in order to show their rela- 
tive impact on the helicity amplitudes, we plot H- = H l 3/HY and Hq = H^/Hq . 
The two approaches agree on the relative smallness of these terms, of a few per- 
cent, across most part of the low q 2 region. Remarkably, this continues to be true 
in the upper part, despite the fact we are neglecting the excited vector mesons 
p*, u* and 0* that populate the region q 2 > 2 GeV 2 . In this regard, we point 
out that these states have a sub-dominant effect in the hadronic structure of the 
electromagnetic current as it is, indeed, probed in e + e~ annihilation experiments. 
Therefore, one would expect their contribution to produce some oscillations in 
the hadronic determination converging, approximately, to the QCDF result at 
larger q 2 . In the numeric computations we ignore these higher-mass resonance 
contributions and we enforce a continuous matching onto the QCDF result at the 
threshold q 2 ~ 2 GeV 2 by means of a suitable g 2 -dependent smearing function 
(see Appendix [B]) . 

On the other hand, in the regions around the resonance poles, the discrepancy 
between the two approaches is maximal and the hadronic model gives contribu- 
tions which are comparable in magnitude to the leading ones. In order to discuss 
this in more detail, consider the schematic form of the semileptonic decay ampli- 
tude (see Appendix pL 



M T = M + e 87C Qvfffv M > : 

(q 2 — m v + im v T v ) 

given by the interference of the leading piece M. including the effects of the 
electromagnetic and semileptonic penguin coefficients 6*7,9.10 and Ai' , the B — > 
VK* amplitude normalized such that \M.\ ~ and e ~ A^ KM ~ 0.05. Then, 

\M T ? = \M\ 2 + e 2 .^ Q l( 2 K :% 2 \Mf 

(q 2 - m v Y + rriyTy 

+ { {q2 ~ m2v)Re[M * M ' ] + mv Tv i ^ m * m '} } ■ 

The contribution of the resonances to the decay rate is sizable and potentially of 
the same order as the leading one only within a g 2 -region of width ~ myTy around 
the position of their poles, in which e x (fvfx* /mv^v) ~ 1- This means, at the 
same time, that integrating over a large enough region of q 2 , Aq 2 ^> myYy, 
suppresses their relative contributions by a factor my Ty / 'Aq 2 . Therefore, we 
conclude that although the effects of the light resonances could alter the line 
shape of observables, they are largely washed out by binning in q 2 . We will see 
this explicitly in the results presented in the next section. 
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3.3 Summary and phenomenological implementation 

In the preceding two subsections we have studied hadronic effects in the contri- 
butions of the various parts of the weak AB = 1 effective hamiltonian to the 
helicity amplitudes. Our main outcome is a strong suppression of the helicity 
amplitudes Hy and H\, which in turn suppresses the coefficients ^3 and Ig in 
the angular distribution. 

• At the factorizable level, this comes about through a double suppression of 
T + by q 2 /m 2 B and A/m&, and a single suppression of V+ by A/m&. In the 
absence of right-handed currents (primed Wilson coefficients) this translates 
to suppressed Hy and H\. 



The contributions of the hadronic weak Hamiltonian enter only into H, 



v 



via h + defined in (23). They are calculable at the leading power in an 



expansion in A/m& in QCDF. At this order, h + = 0. 

• Long-distance "charm loop" effects can be studied in an expansion A 2 / (4m 2 — 
q 2 ). At zeroth order, one recovers the non-factorizable form factor QCDF 
expressions, which give vanishing contributions to h + . The first nonzero, 
long-distance, contribution to h+ arises at order A 3 / (4mbm 2 ) relative to the 
amplitudes Hy and H v . Long-distance contributions of the chromomag- 
netic operator Q 8g are expected to be small, with the contribution to h + 
suppressed (at least) as A 2 /m 2 . 

• The remaining effects are suppressed by small CKM and Wilson coefficients. 
We estimate long-distance corrections to QCDF by means of a resonance 
(VMD) model, which embodies a suppression of the contribution to h + , as 
a consequence of the helicity structure in B — > VK* decay. After binning 
the model predicts very small corrections in the other amplitudes h- and 
h . 

In our phenomenological analysis, we employ our model-independent parameter- 



isation of form factors through (55) and (56) (keeping terms up to and including 
0(q 2 /m 2 B )), together with the numerical values and ranges given in Sect. 3.1 For 



the contributions of the hadronic weak Hamiltonian, we employ the leading-power 



QCD factorization expressions 23,29 . To estimate the long-distance corrections, 



we employ the following model for long-distance charm loop effects, 

/i_| cg , LD = O.le^Cf 1 , (76) 
^+Uld = 0.02e^+C 7 SM , (77) 

where <fi± are arbitrary soft rescattering phases, comprising a conservative inter- 



pretation of the numerical findings of 46 and the hierarchy h + h_. We will 



also allow for an extra long-distance contribution /io|cc,ld = 0.2e^°Cf M , as there 



30 



is no power suppression of ho] this should be considered an ad hoc model but 
does not impact on the observables emerging as clean in the phenomenological 
analysis below, which only involve helicities ±1. We have increased the magni- 
tude of these effects beyond the error estimates of 46 , such as to accomodate 



within it the small long-distance contributions from the chromomagnetic penguin 
operators. For the light-quark Hamiltonian, we estimate possible long-distance 
corrections by means of the model described in Sect. 3.2.3 Input parameters are 
summarized in Appendix |A. 2 



4 Phenomenology at low q" 

4.1 Observables for B K*i + i~ 



In general, there are twelve g 2 -dependent observables (shown in 39 neglecting 
tensor operators, but this remains true in the their presence) that are accessible 
through a full-angular analysis of the B — > K*£ + £~ decay rate and which corre- 



spond to the angular coefficients h(q 2 ) in Eq. (38). In the absence of scalar and 
tensor operators, which includes the SM, Jg = 0, and there is one relationship 
among the remaining coefficients, reducing the number of independent observ- 



ables to ten 45 . If one furthermore assumes rri£ = 0, two more relations can be 
established, 

3/is = hs, he = —he, (78) 

leading to eight independent observables. 

The analysis of the CP-partner decay B — > K*£ + £~ gives a same amount of 
independent observables as in the B decay, the Jj's. In this sense, it is useful to 
define the following combinations of Jj's and Jj's, 

£< = A, = (79) 



which can be used to construct a variety of CP-averages and asymmetries 37,39 . 



4.2 The B° K*°fi + f£~ decay 

In Fig. [6j we show the SM predictions for the eleven angular coefficients available 
in this case and normalized by the B° decay rate at low q 2 . The solid (red) and 
the dashed (green) lines correspond to the prediction including the light-quark 
contributions in the hadronic model or in QCD factorization, respectively (see 
Sec. [3]). The inner (red) error band is the uncertainty derived from the hadronic 
parameters (soft form factors, decay constants,. . .), the CKM parameters and the 
renormalization scale. The intermediate (blue) and outer (green) bands result 
from the addition in quadratures of the unknown factorizable and charm-loop 
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Figure 6: Angular coefficients normalized by the B° decay rate around the 
low-g 2 end-point. Solid (red) and dashed (green) lines correspond to the SM pre- 
diction including the light-quark contributions in the hadronic model or QCD fac- 
torization, respectively. The error bands stem from the hadronic and CKM uncer- 
tainties and renormalization scale dependence (inner-red), plus the unknown fac- 
torizable power-corrections (internal-blue) and plus the non-factorizable charm- 
loop uncertainty (outer-green). The errors are added, subsequently, in quadra- 
tures. 
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Figure 7: Differential branching fraction, and the "clean" observables 
around the low-g 2 end-point. We show in black the experimental results for the 
two first observables in the bins [0.05, 2] GeV 2 and [2, 4.3] GeV 2 (6j. The color 
code is as in Fig. |6j 



power corrections, subsequently. Factorizable corrections are estimated using 
Eq. (56), and the charm- loop uncertainty are modelled according to Eqs. (77). 

The main source of uncertainties in the Ij's stem from the soft form factors 
and, in some cases, from the charm-loop. In particular, for the coefficients pro- 
portional to Hy A , J 3 and Ig, the latter source is, by far, the most important. On 
the other hand, it is remarkable that the uncertainties in the coefficients arising 
from the unknown factorizable power corrections are negligible at low q 2 . This 
effect is due to the constraints imposed by the exact relations (44). Finally, notice 
that the vector-meson resonances alter significantly the line shape of most of the 
Jj's, except for those oc Hy A due to the suppression of the corresponding helicity 
amplitude in the B — > K*V decays (see Sec. K3I). 
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4.2.1 CP-averages: The branching fraction, Fl and the P-basis 

Each of the observables that can be constructed out of the CP combinations in 



Eqs. (79) has a different sensitivity to the various standard and non-standard 



Wilson coefficients. In order to maximize these sensitivities, it is important to 
find a set observables with reduced dependence to the uncertain hadronic param- 
eters underpinning the theoretical predictions, in particular the B — > K* form 
factors 32 , 38 . Indeed, a proper set of "clean" observables can be obtained using 
suitable ratios of angular coefficients 53 57 66 . Following this strategy, we use 



the following set of CP-averaged observables 66 



Pi 



2S 



P? 



2.s 



8S 



2.s 



4S 



2.s 



P' 



p' 



P' 



-Y>2s^>2c 2y— S2 S S2c 2y— S2 S S2c 

These observables remain independent of each other when the leptons are as- 
sumed to be massless. As it was explained above, two more independent observ- 
ables can be defined in case one lifts this assumption and incorporates lepton-mass 
dependent contributions to the decay rate. These other observables have some 
special features and we discuss them separately in Sec |4.2.3| Nevertheless, we 



keep all the lepton-masses at their physical values in the differential decay rate 
and, in particular, in the expressions of the angular coefficients. 



The set of observables in Eqs. (85) is considered clean in the sense that the 



soft form factors cancel at LO in the ratios so that the uncertainty stemming 
from the form factors is suppressed in a s (as computed at NLO in QCD factor- 
ization) or A/rrib (unknown factorizable power corrections). Notice that some 
of these observables can be related with others previously defined in the litera- 

in 



ture 32,38,53 . In order to form a complete set of independent observables, 



the (pseudo) scalar-less case, we add the decay rate V and Fl, 



V 



ldT + dT 
2 



dq 2 



lc 



->2c 



2(3S ls -E 2s )) 



(80) 



3S ls 



J 2s 



3S lc 



J 2c 



satisfying Ft — 1 



2 p 4 p v j 

■Fl- Another option would involve, e.g., the forward-backward 
asymmetry (41) although notice that it can be obtained straightforwardly com- 
bining P2, Fl and the decay rate. 

In Fig. [7l we plot the differential branching fraction, Fl and the clean ob- 
servables P{_ 6 , and we follow the same color and line code as the one used in 
Fig. |6j Also, we show in black boxes the experimental results for the two first 
observables in the bins [0.05, 2] GeV 2 and [2, 4.3] GeV 2 that have been measured 
by the LHCb collaboration (6l. These measurements agree very well with the SM 
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predictions. By comparing the first row with the second and third rows of the 
panel in Fig. [7], we ratify the advantage of using a set of observables with reduced 
theoretical uncertainties [32j[57j[66] . While for the differential branching fraction 
and Fl the limited knowledge of the hadronic parameters (specially the soft form 
factors) is the dominant source of uncertainty, for the Pj the power corrections 
become much more important. In the latter case, the enforcement of the form 



factor relations (44) is essential to constrain the size of the factorizable power 
corrections at very low q 2 , cf. Pi and P3. As a consequence, the uncertainty 
in this region is dominated by the power corrections to the charm loop. This 
has to be interpreted as an outcome of the not-very-precise knowledge we have 
of these contributions, as compared with the information we can gather for the 
form factors. As for the effect of long-distance effects in the light-quark contri- 
butions we see that, indeed, they can modify abruptly the line shape of most of 
the observables in the neighborhood of the vector-resonance poles. Nevertheless, 
these effects are again tiny in Pi and P 3 (which are oc Hy A ) and, as argued in 
Sec. [3j they dilute after binning in q 2 . Moreover, Pi and P3 as pure helicity-1 
objects are free from S-w&ve contamination. They emerge as especially clean null 
tests of the Standard Model. 

On the experimental side [6], what is measured are ratios of binned observables 
(rather than binned ratios), hence it becomes necessary to define the observables 
not as g 2 -integrals of functions of the Jj's but rather as the same functions of the 
corresponding binned angular coefficients; only in case the coefficients are slowly 
varying functions of q 2 in the bin considered, the two methods give approximately 
the same result. As one can deduce from Figs. [7] and [6j this is not the case around 
the low-g 2 end-point. With these considerations, we define the following set of 



CP-averaged binned observables 66 



2 \ dq 2 I ' v ' 4(F) 



2(E 2s )' v Zl 8<E 2s )' v " 4<E 2s )' 



v/-(S 2s )(S 2c )' 2v/-(E 2s )(E 2c )' 2 v /-(E 2s )(E 2c )' 

(82) 



where (E;) = j q h ax Zi(q 2 )dq 2 . In Table 



we present the predictions for the 



binned observables, as defined in Eq. (82), in some bins of interest 111 . In 
particular, we highlight the results on the low-g 2 bin [0.1, 1] GeV 2 , which presents 
an overall small theoretical (relative) error, as compared with the higher q 2 bins, 
due to the small form factor uncertainties in this region. It is also remarkable that 
the estimate on the contributions from the light hadronic resonances comes to 



35 



Table 3: Results and error budget on the binned CP-averaged observables of the 
muonic mode. 



Obs. 



10 7 x (f ; 



\a 2 ■ a 2 1 

[0.1, r 

[0.1, 2 
[2, 4.3 
[1,6] 



Result 



-0.23 
-0.20 
-0.39 
-0.38 
-0.33 
-0.26 
i c+0.8 

L -°-o.e 



o.8i; 
1.13: 

0.62: 



Hadronic Fact, c-quark Light-quark 



+0.20 
-0.17 
+0.36 
-0.24 
+0.27 
-0.21 
+0.6 
-0.5 



+0.03 
-0.03 
+0.08 
-0.07 
+0.19 
-0.15 
+0.46 
-0.37 



+0.10 
-0.10 
+0.13 
-0.12 
+0.02 
-0.01 
+0.05 
-0.05 



±0.00 
±0.02 
±0.00 
±0.02 



(Fl) 



[0.1, 1 
[0.1, 2 
[2, 4.3 
[1,6] 



0.20^ 



-0.11 
-0.10 

u - oi -0.12 

0.7018:11 



+0.10 
-0.09 
+0.15 
-0.11 
+0.09 
-0.13 
+0.11 
-0.13 



+0.02 

-0.02 

+0.04 

-0.04 

+0.07 

-0.9 

+0.09 

-0.11 



+0.03 
-0.02 
+0.04 
-0.03 
+0.02 
-0.02 
+0.02 
-0.02 



±0.01 
±0.01 
±0.00 
±0.00 



10 2 x (P x ) 



[0.1, 1 
[0.1, 2 
[2, 4.3 
[1,6] 



9 q+ 3 - 2 

z - y -3.1 
o n+3.5 
O.U_ 3 .4 
_1+7 
L -5 
o+8 
z -6 



+0.8 
-0.1 
+0.8 
-0.2 
+1.6 
-0.8 
+1.3 
-0.8 



+1.2 
-1.3 
+1.7 
-1.7 
+7 
-5 
+8 



+2.9 
-2.8 
+2.9 
-2.9 
+1.8 
-1.6 
+1.6 
-1.4 



±0.0 
±0.1 
±0.0 
±0.0 



10 x (P 2 ) 



10 2 x (P 3 ) 



[0.1, 1 
[0.1, 2 
[2, 4.3 

[1,6] 
[0.1, 1 
[0.1, 2 
[2, 4.3 

[1,6] 



+0.15 
0.17 
0.19 
0.26 

1+1.4 
°- i -1.6 

1 4 +1.5 
-'-•^-1.5 

1.5 
1.2 
6 
3 
1.2 
1.2 

n 3 +1 -° 

u -°-1.0 



1.02 
1.57 



-0.1 
-0.21} 
-0.3 



+0.08 
-0.13 
+0.08 
-0.20 
+0.8 
-0.8 
+0.8 
-0.7 

+0.0 
-0.2 
+0.0 
-0.2 
+0.1 
-0.3 
+0.1 
-0.3 



+0.10 
-0.09 
+0.13 
-0.13 
+1.0 
-1.2 
+1.2 
-1.1 

+0.1 
-0.1 
+0.1 
-0.1 
+0.7 
-0.8 
+0.6 
-0.6 



+0.08 
-0.07 
+0.11 
-0.10 
+0.5 
-0.7 
+0.5 
-0.6 

+1.5 
-1.2 
+1.6 
-1.2 
+1.0 
-0.9 
+0.8 
-0.7 



±0.00 
±0.04 
±0.0 
±0.0 

±0.0 
±0.0 
±0.0 
±0.0 



10 x {Pi) 



[0.1, 1 
[0.1, 2 
[2, 4.3 
[1,6] 



-5.1" 



-0.9 
-0.4 
o 0+1.4 
-O.8_o.6 

4.61^ 



+0.8 
-0.0 
+1.2 
-0.0 
+0.9 
-1.4 
+0.9 
-1.2 



+0.2 
-0.2 
+0.3 
-0.3 
+1.3 
-1.4 
+1.1 
-1.2 



-0.3 

-0.2 

-0.4 

-0.4 

-1. 

-0.9 

-0.8 

-0.8 



±0.3 
±0.4 

±0.0 
±0.0 



10 x <i*> 



[0.1, 1 
[0.1, 2 
[2, 4.3 
[1,6] 



6 q +0 - 8 

u - y -0.5 



;+0.7 
-1.1 
-3.1 
-2.7 
o o+3.0 
-^•S-2.6 



5.5: 
-2.5: 



+0.6 
-0.1 
+0.1 
-0.8 
+1.5 
-1.1 
+1.3 
-1.1 



+0.3 
-0.3 
+0.6 
-0.6 
+2.5 
-2.3 
+2.5 
-2.2 



+0.2 
-0.2 
+0.3 
-0.3 
+0.9 
-1.0 
+0.8 
-0.9 



±0.3 
±0.2 
±0.0 
±0.0 



10 x (P^) 



[0.1, 1 
[0.1, 2 
[2, 4.3 
[1,6] 



0.7 
0.8 
0.7 
0.8 

0.91 1 - 



-0.8 
-0.8 



-1.0 

-0.71?:* 



+0.3 
-0.4 
+0.2 
-0.5 
+0.3 
-0.5 
+0.2 
-0.5 



+0.0 
-0.0 
+0.0 
-0.0 
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-0.1 
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-0.1 



+0.7 
-0.7 
+0.7 
-0.7 
+0.9 
-0.9 
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±0.0 
±0.0 
±0.0 
±0.0 
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Figure 8: Differential Aqp and P^ p at low q 2 . For Acp we show the experimental 
data point at the low-g 2 bin [0.05, 2] GeV 2 . The charm-loop uncertainty and 
factorizable power corrections are negligible in this case. 



be, for the CP averages, negligible. In this regard, the largest pollution appears 
in the bins including the 0(1020) resonance and in the observable P' A . As we will 



see in Sec. 4.4 low-g bins are also important for their sensitivity to the Wilson 



coefficients C 7 and C' 7 . 

4.2.2 CP asymmetries: A C p and P 3 CP 

The CP asymmetries are observables suitable for searching for sources of CP- 
violation beyond the CKM mechanism. In B — > K*£ + £~, a weak phase arises 
from the interference between the contributions weighed by V t bV t * s and V u bV* s , 
so all the CP-asymmetries for this decay in the SM are suppressed by a factor 



\ ut = V ub V* s /V t bV t *, which imaginary part is of order ~ fj\ 2 ~ 10~ 2 18 . Another 
important observation is that the sensitivity to a CP- violating phase is modulated 
by a factor sin 5 S sin 5 W or cos 5 S sin 5 W depending on whether the asymmetry is 
T-even or T-odd (i.e odd under the transformation <p — > —<p). In this sense, QCD 



factorization predicts the strong phases, 5 S , to be very small 23 29 and, conse- 



quently, the latter asymmetries have been specifically singled out as interesting 



observables for the detection of new CP-violating phases beyond the SM 37 . 

Nonetheless, a hadronic, rather than a partonic treatment, is suited for the 
description of the contribution of the light-quarks at low-g 2 (see Sec. [3]). In the 
context of VMD that is developed in this work, we have concluded that the 
overall contribution to the CP-averaged observables is negligible and consistent 
with the intrinsic suppression induced by small Wilson coefficients or by the CKM 
factor X ut . For the CP-asymmetries, though, large differences between the two 
approaches can be expected since the leading contribution is given by the latter 
pieces. To be more precise, large strong phases can be obtained in the hadronic 
picture as they naturally arise from rescattering and "dressing" of the resonance 
poles. Therefore, at low q 2 a suppression oc sin 5 S ~ of the T-even asymmetries 



37 



cannot be expected and sensitivity to the weak phase (s) in these observables can 
be generated by long-distance effects. Conversely, a factor cos 5 S close to 1 is not 
guaranteed on general grounds at low q 2 , which can hinder the efficient use of the 
T-odd observables for probing CP-violation in this regime^] 

As an illustration, we show in the left-hand side of Fig. [8] the CP-asymmetry 

. T-f 3A lc - A 2c + 2 (3A ls - A 2s ) 

Acv = fTf = 4F ' (83) 

which is T-even and, on the right hand-side, the T-odd asymmetry P^ p |l!2| , 

T = ~ (84) 
for the decays of the neutral mesons in the muonic mode. The latter observable is 



a redefinition of A 9 18,39 following the same "cleanness" principle applied to set 



the P-basis for the CP-averaged observables 112 . Aqp has been measured at the 



.B-factories [4 10 and, recently, binned results and the most precise measurement 
to date of the total have been provided by the LHCb collaboration |113| . An 
interesting outcome of the latter measurement is that the total (integrated) Acp 
is negative and of order ~ 5%, with a sizable and negative contribution stemming 
from the lowest-g 2 bin. This seems to be in contradiction, at la level, with the SM 
prediction as it is obtained in QCD factorization. On the other hand, this type 
of substantial CP-asymmetries could be the smoking gun for some NPs scenarios 
(see e.g. (54]). However, and as shown in the left-hand side panel of Fig. [8] long- 
distance contributions of the light quarks can have an important effect in Aqp- 
In fact, we obtain that the uncertainty produced by these can be rather large, 
cf. Aqp = (0^)% in the [0.05, 2] GeV 2 bin, with an error band one order of 
magnitude larger than in QCD factorization. 

As for P3 p we observe that it is practically zero in the SM. This observable 
(in the form of A$) has been tagged as a benchmark for the detection of a non- 



standard weak phase that might surface through the Wilson coefficient C 7 37 39 
and our findings show that it remains a clean null test, safe from long-distance 
charm contributions and contamination by light resonances as a consequence of 
being oc Hy A . We investigate the sensitivity of this observable, together with 
other combinations of the angular coefficients ^3 and I3, to BSM values of C 7 , in 
Sec. HI 



In any case, we want to stress that a full treatment of the hadronic uncertainties on the 
CP-asymmetries is beyond the scope of the present work. For instance, neglecting the effect of 
higher-mass p, to and 4> resonances present across all the lowest-g 2 region might not be as safe 
approximation for the CP-asymmetries as it is for the averages. 
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4.2.3 Non-vanishing lepton masses: The observables M\ and M 2 

Two more independent "clean" observables can be denned for the muonic mode 
at low q 2 [57], Q 

M > - 1 Ph. - (2 + P)h s , 1 P 2 h c + h c 

1_ l-/3 2 4/ 2s ' M *~ \-p I 2c ' (85) 

which involve the pieces proportional to 4m 2 jq 2 present in Ii s and I\ c , 
M[ -- 





2 + 




2 - {V -> A)) 




2 + 




2 + (V -> A)) 



g 2 /(2m 2 )(|i/p 


2 + /3 2 


|#S 


2 ) + (\Hy 


2 _ 




2 ) 


(1^1 


2 + 




2 ) 



^ = ^ ^ M ; ' ~ ' "'^ (86) 



The two questions that naturally arise are, first, what can we learn from a 
measurement of these observables; second, how these pieces enter in the formulas 
of the differential decay rate expressed in terms of a given basis of observables, 
which are usually defined in the me — > limit. To start off, notice that the 
contribution of these pieces to the differential decay rate vanishes except in the 
vicinity of the low-g 2 end-point, e.g. 4m 2 /(0.5GeV 2 ) ~ 0.09. In this region it 
then follows that 

M[ = \ + 0(q 2 /m 2 B ) 2 , M' 2 = l + 0{q 2 /m 2 B ). (87) 

To obtain the first relation it is sufficient to note that the photon pole in Hy 
makes to be of 0(q 2 /m 2 3 ) in comparison. Thus, this relation is completely 
model-independent, in the sense that it is insensitive to the long-range hadronic 
uncertainties as well as to the particular short-range structure of the decay. In 
order to derive the second relation, we first omit the effect of (pseudo) scalar 
operators. Then, using the corresponding expressions for the amplitudes it is 
straightforward to obtain 

\H° V \ 2 + q 2 /(2ml)\H P \ 2 - \H° A \ 2 = \H° V \ 2 + \H° A \ 2 + Rem., (88) 



where 

N 2 X 
2m 2 B q 



Rem. = ^^|C 10 - C[ \ 2 (S(q 2 ) - V (q 2 )) ~ 0(q 2 /m 2 B ) (89) 



due to the exact form-factor relation (45). The relation for M' 2 is not model- 
independent as it can be broken by contributions from BSM scalar and pseu- 
doscalar operators, although the size of these effects is very constrained by recent 



9 For the sake of making our discussion in this section more transparent, we slightly redefine 
these observables as introduced in that reference, factorizing out a Am 2 /q 2 piece. 
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Figure 9: Results for M[ 2 and Ri^. The color code is as in Fig. [6j whereas the 
black line is the suppression factor 4m 2 /g 2 . 



measurements of the B s — > n + fi~ decay rate. Indeed, if one uses the experimen- 
tal upper bound given by the LHCb collaboration at 1-a level [I] and the SM 
prediction in Ref. [2] to constrain the size of these pieces (assuming C10 = Cfo 1 ), 
one obtains that the relation for M' 2 is broken only at the ~ 1% level. 



The Eqs. (87) mean that, in practice, M[ 2 are not independent observables 
in the q 2 region in which the decay distribution has sensitivity to their effects. 



Taking this into account, one can generalize Eqs. (78) with expressions that are 
valid at finite me and then, also, close to the low-g 2 end-point. Namely, we define, 

ha = ^(4 - (3 2 )hs, he = ~{2 - /3 2 )/ 2c , (90) 



that convert into the conventional formulas (78) in the massless limit. In the 
upper panel of Fig. [9] we show the q 2 dependence of M[ and M' 2 at low q 2 using 
the same color code as in Figs. [6] and [7} We superimpose in black the value 
of the factor 4m M /g 2 suppressing the contribution of these pieces to the decay 



rate. In the lower panel we explore the validity Eqs. (90) by showing the ratios 
Ri,2 — hs,2s/hs,2s- As one can check in the plots, these relations work at better 
than 1% of accuracy in the region of interest. 



4.3 The B u -» iT°e + e decay 



In comparison with the muonic mode, the electronic decay is interesting be- 
cause its low-g 2 end-point is at q 2 = 4m 2 ~ 10~ 6 GeV 2 . Hence, it presents an 
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Figure 10: Differential branching fraction and Pi for the electronic mode at low 
q 2 . The color code is as in Fig. |6j 



Table 4: Results and error budget of the integrated branching fraction and Pi 
observable in the [0.0009, 1] GeV 2 bin for the electronic mode. 



Obs. 


Result 


Hadronic 


Fact. 


c-quark 


Light-quark 


10 7 x (f ) 


9 A q+0.66 
z -^°-0.47 


+0.50 
-0.39 


+0.10 
-0.05 


+0.42 
-0.25 


±0.03 


10 2 x (Pi) 


o 7+3.0 
Z. 1 _ 2 7 


+0.8 
-0.1 


+1.0 
-1.2 


+2.7 
-2.3 


±0.0 



enhanced sensitivity to the physics associated to the photon pole becoming a 



golden mode 
operator Q$ 
tronic mode 



xxr probing possible BSM effects related to the magnetic penguins 
. On the experimental side, B— factories have observed the elec- 
and prospects at LHCb are discussed in |114| . This will also be 
an interesting target for super flavour factories. Our detailed error assessment of 
the low-g 2 region in the muonic mode revealed that the low-g 2 end-point is actu- 
ally the least prone to form factors uncertainties due to the constraints (44), at 
the same time as the hadronic contributions, for CP-averaged observables, come 
out to be of negligible size. Consequently, the B° — >■ K*°e + e~ decay emerges as 
highly relevant. 

In Fig. 10 we show the SM predictions for the differential branching fraction 
and the Pi observable, which is especially sensitive to NPs effects in the magnetic 

In Table |4| we show the integrated results for the bin 

As one can see, the integrated 
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penguin operators 32 , 57 
[0.0009, 1] GeV 2 , that is the one proposed in 
branching fraction is larger than the one for the muonic case in its lowest bin. 
Interestingly enough, our result agrees with the estimate B ~ 2.2 x 10 that is 



obtained assuming that the decay is entirely driven by the photon pole 21,114 



This dominance of the transverse amplitudes in this integrated decay rate is also 
reflected by a very small integrated longitudinal polarization in the same bin, 
namely (F L ) = 0.077±g : gg. 
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Figure 11: Study of the sensitivity of the observables Pi and P% p to the a purely 
real or purely imaginary NPs contribution through C 7 (dashed line, blue bands). 
These are confronted with the SM expectation (solid line, red band). 



4.4 Sensitivity to C' 7 

The analysis of the low q 2 region of the B — > K*£ + £~ can provide tight constraints 
on NPs scenarios with right-handed flavour-changing neutral currents, specially 
those giving contributions to the chirally-flipped magnetic penguin operator 0' 7 . 
This is due to the fact that the angular coefficients ^3 and Ig, at low-g 2 , are 

I 3 oc Re (HyHf) , Ig (x Im (HyH+*) , (91) 

where Hy oc C' 7 /q 2 , so, approximately, they vanish unless C' 7 7^ 0. (Corrections 
involving H\ are also suppressed by the smallness of C' g and C[ in the Standard 
Model, but any BSM effects generating H\ are suppressed at q 2 ~ due to the 
absence of a photon pole.) In the SM, small contributions to these observables 
are generated by the strange-quark mass and other effects quantified in this work 
as contributions to the Hy helicity amplitude. Other decays and observables 
provide valuable and independent constraints on the C 7 and the C' 7 planes, in 
particular the inclusive B — > X s 7 decay and the isospin and the time-dependent 
CP-asymmetries in the exclusive B — > K*^ decay (see e.g. 64|). The interest 
of the radiative decays onto higher-mass K* resonances has been also recently 
pointed out [48 1. In this work we focus on studying the sensitivity of the vicinity 
of low-g 2 end-point to the chirally-flipped Wilson coefficient C 7 . A more compre- 



hensive analysis should also consider studying new-physics effects in C7 63 -65 
although ^3 and Ig can be used to efficiently constrain also this Wilson coefficient 
only if C 7 is far from zero. 

In the context of CP-combinations, one can construct 4 independent observ- 
ables with these angular coefficients and their CP-conjugates. However, 7 3 and Ig 
are a CP-odd and a CP-even observable, respectively, and the combinations A3 
and Eg become not very sensitive to the chirality of dilepton pair. Therefore, only 
1 CP-average and 1 CP-asymmetry, constructed from E3 and Ag, in order, are 
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sensitive to either the real or imaginary parts of C 7 '\ The corresponding observ- 
ables in the P-basis, P\ and P% p , are a convenient choice for their intrinsically 
reduced theoretical uncertainty. In order to investigate the sensitivity of these 
observables to C' 7 , we investigate two a priori fictitious scenarios: the first one in 
which C' 7 is real and with a magnitude of 10% the value of C 7 {rrib) in the SM, 
and second where C 7 is a pure imaginary number with and absolute value of 1% 
the value of C^lnib) in the SM. The result for Pi in the first scenario is plotted 
as a dashed line confronted against the SM prediction represented as a solid line 



on the left panel of 11 The corresponding bands are the overall uncertainties in 
these scenarios calculated as those presented for the SM in the previous sections. 
On the right-hand side, we plot P^ p for the second scenario compared, again, 
with the SM prediction and following the same line and color code as on the left 
panel. 

The first remarkable outcome of this analysis is that, with the current theoret- 
ical uncertainties, Pi is sensitive to a C 7 with a real part of about 10% x C 7 M (rrib) 
only in the region between q 2 = 4mj and q 2 ~ 3 GeV 2 . The same conclusions 
apply to P 3 CP although, in this case, the attainable accuracy is much higher as it 
can be sensitive to complex phases with an absolute value below 1% x C 7 M (rrib)- 
This, in principle, is an important result as it singles out P|" p as a theoretically 
extremely clean observable to constrain (or favor) BSMs generating new CP- 
violating phases through C 7 . This occurs despite the cautionary remarks made 
in Sec. 4.2.2 about the CP-asymmetries. In fact important interference with 
strong phases, with a drastic reduction of the sensitivity, seems to appear only in 
the immediate surroundings of the resonance poles. Notice, however, that P^ p 
is an special observable as it vanishes in the SM so the cleanness of other CP- 
observables can not be concluded without a devoted study on their corresponding 
hadronic uncertainties. 



5 Conclusions and outlook 

We have performed a comprehensive analysis of angular observables in the decays 
B — > K*£ + £~, £ = fi, e, paying particular attention to non-factorizable hadronic 
uncertainties. We exploit a suppression of the positive-helicity amplitudes in the 
Standard Model that, as we have shown, still holds when taking into account 
long-distance charm and light-quark non-factorizable effects in a conservative 
way. As a result, we established that of those observables that are "clean" in 
the factorizable approximation in the large-energy/heavy-quark limit, the two 
observables Pi = Aj) and P^ p are found to nearly vanish in the SM with very 
small errors in the presence of the long-distance non-factorizable corrections. At 
the same time they are highly sensitive to right-handed currents. Importantly, 
we find that the lower end of the low-g 2 range is not only theoretically clean, but 
provides the best sensitivity to the Wilson coefficient C 7 , down to a (theoretical) 
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limit of below 10% for the real part, and as low as 1% for the imaginary part. This 
also raises the profile of the decay into electrons, I = e, which is dominated by the 
photon pole, due to the small electron mass. We feel that this process deserves 
further attention. For other observables, which do not vanish for Hy = 0, the 
theoretical control is more compromised by long-distance effects. Improving it 
would seem to require further theoretical advances. 

As a further result, we have found a way to generalise the two well-known 
relations 3Ii s = I2 S and I\ c = —li c , which hold for vanishing lepton mass, to 
the massive case. This may be of use in the experimental determination of the 
angular distribution. 

The formulation in terms of helicity amplitudes allowed us to show, by adapt- 
ing the LCSR formalism of (46], a strong A 3 /(4m 2 m;,) power suppression of the 
long-distance contribution /i+| c g,LD from the charm part of the hadronic weak 
Hamiltonian entering in the positive-helicity amplitude Hy, which vanishes at 
leading power in QCDF. Analogous contributions from the chromomagnetic op- 
erator are suppressed by A 2 /m 2 . Likewise, we exclude the possibility of non- 
negligible, light-resonance long-distance contributions to Hy once the (A/m^) 2 
suppression of the corresponding B — > K*V helicity amplitude is taken into 
account. More generally, we observe that, once data are binned, these light- 
resonance contributions are small in all observables, even though this is not true 
of the unbinned quantities. On the other hand, the amplitudes Hy and H v can 
receive important long-distance charm-loop contributions, at a level of (9(10%) of 
the leading-power QCDF result, as is clear from the numerical values in [46|. We 
have also argued that the much smaller estimate of these contributions of |34| is 
likely to suffer very large, incalculable corrections. This limits the theoretical con- 
trol over the other "clean" observables, beyond P\, P3, and their CP asymmetry 
counterparts, that, although independent of H v , involve Hy. 

Our phenomenological analysis makes certain assumptions about the values 
of helicity form factors and nonfactorizable corrections to the helicity amplitudes. 
To get a better quantitative understanding of the theoretical errors, it would be 
desirable to have direct calculations of form factors in the helicity basis. This 
should be straightforward within LCSR, or other methods. It is also important 
to stress that we have not derived a sum rule for /i+| C c,ld, but only shown a 
parametric suppression of it (and of the Q$ g long-distance contributions.) An 
attractive aspect of the LCSR method is that it avoids endpoint divergences al- 
together, and a sum rule allowing for a quantitative estimate could in principle be 
obtained. We feel that that is a path worthwile pursuing. On the other hand, for 
the light-quark non-factorizable effects it would be interesting to go beyond the 
resonance model in assessing the size (and helicity structure) of possible correc- 
tions to the QCDF results, perhaps by combining QCDF, LCSR, and dispersion 



relations along the lines of 94 . While this is unlikely to have a strong impact on 
P\ or P 3 CP , it may be more important in other observables, and in particular for 
the direct CP asymmetry Aqp, which, as our resonance model suggests, could be 
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relatively sensitive to any large, "soft" strong phases associated with light-hadron 
intermediate states. 
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Note added. After submission of this work to the arXiv, a paper including 
a new LCSR approach to the matrix elements of Qg g appeared 122 . It would 



be interesting to adapt this method to the approach of the present paper. We 
also note work on the effect of tensor operators in B — > K*£ + £~ in the high-g 2 
region 
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A Notations and conventions 



A.l Polarisation vectors 

For a particle moving in +z direction with four-momentum k 11 = (E, 0, 0, |k|) we 
use 



e"(±) 
e"(0) 



1 

i 



(0,Tl,-«,0), 
(|k|,0,0,£). 



(92) 
(93) 
(94) 



In parentheses is the helicity, we sometimes alternatively put this in a subscript. 
The subscript t denotes "timelike" polarisation, which is a spin singlet. The sign 
conventions ensure that the helicities ±, form a triplet obeying conventional 
relative signs and normalisations, e.g. [J + ] M jy e t/ (0) = v2e M (l) . 

The polarisation vectors satisfy the normalisation and completeness relations 



6 = ±,0, 



c a b,fi 



9ab, 



9ab 



I- 

0. 



a = 
a = b = t. 
otherwise 



(95) 
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A=±,0 o,6=t,±,0 

For any other direction n, we define the polarization vectors by applying a stan- 
dard rotation after boosting to the "rest frame" where = 0,0,0). In 
that frame, 



^{h- A) 



e -i<pJz e -i0J y e +ic/>Jz 



e»(z; A) = J] 4a We i(A ~ V) A'), (97) 



where 0,0 are the spherical coordinates of n and d 3 m , m (9) is the matrix element 
of the rotation operator e~ l6Jy (conventions as in 115| ). For arbitrary direction, 
the final form is usually more convenient. 

For h = —z (ie 9 = ir) we rotate at <ft — 0. This means 

e(-z;±) = e(z;=F), e{-z; 0) = -e{z; 0). 

Boosting back to the laboratory frame, we have 



e"(-*;0) 



v/2 
1 



(£,0,0,-|k|), 



(0,±l,-i,0) = ^(z;T), 
|k|,0,0,-£). 



(98) 
(99) 
(100) 



A. 2 Standard Model parameters 



Table 5: Some of the SM parameters used in this work (in GeV). 



M w m t {m t ) m bjPS m c 


m s {2 GeV) 


A (5) 


80.4 172(2) 4.8(2) 1.3(1) 


0.094(3) 


0.214(8) 



Table 6: CKM parameters in the Wolfenstein parameterization used in this 



work 116 



A 


A 


P 


Tj 


0.22543(8) 


0.805(20) 


0.144(25) 


0.342(16) 



We use the PS-subtracted definition for the 6-quark mass [23] . This is obtained 
from the MS-definition (m&) through the pole mass definition [rrib) using the 
subsequent identities 



rh h {n) 



m b [l+ 



3 log 



- 4 

:>, /'/ 



A' 



+ 0(al 



(101) 
(102) 
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Table 7: Wilson coefficients of the SM at \i 
to NNLL accuracy. 



4.8 GeV, in the basis of 73 and 





c 2 


c 3 


c 4 


c 5 


c 6 


Cf 




c 9 


Cio 


-0.144 


1.060 


0.011 


-0.034 


0.010 


-0.040 


-0.305 


-0.168 


4.24 


-4.312 



The PS-mass at pf = 2 GeV is then identified with the pole mass, mb pg(2 

GeV) = 4.78+^;o? ~ 4.8(2) GeV [TT7| . We use m b as the potential-subtracted 
definition unless it is stated otherwise. For the strange quark mass appearing at 
tree level we use the result obtained from LQCD calculations, rh s {2 GeV) = 94(3) 
MeV 
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Other SM parameters used in the calculation are listed in Tables |5j 
and [6] For a s (p) we use the three-loop renormalization-scale evolution 119 , for 
a e m(Q 2 = 0) ~ 1/137 and for G F = 1.16637(1) x 10~ 5 GeV" 2 . Finally, we show 
in Table 7\ the values of the Wilson coefficients of the SM at p = m& GeV, in the 
basis of 73 1 and to NNLL accuracy. 



A. 3 Leading-twist LCDAs 



We use the parameterization of the light vector-meson LCDAs of Ref. 77 in terms 
of the wave functions &K*,a{u) and the semileptonic decay constants. However, 
we use the conventions and values of Refs. 29 , which have a sign difference in 



the odd-coefficients as compared with those of the LCSR calculation. Namely, 



K*,a 



Quu 



l + 3ai, a £ + a 2 , a -(5£ 2 - 1) 



(103) 



with u = (1 — it), £ = 2m — 1 and a^ a Gegenbauer coefficients. There are two decay 
constants for a light vector meson. The longitudinal constant can be obtained 
from experiment using the branching fractions for the electromagnetic decays 
of the neutral vector resonances and applying SU(3)p symmetry. We obtain 
f K * = 220(5) MeV using the p decay and so this value is consistent with the 
one obtained from sum-rule calculations [77]. Sum rules also predict the tensor 
(transversal) decay constant, /k-*,_l(1 GeV) = 170(20) MeV. The latter is consis- 
tent with the LCSR determination f K *,±( 2 GeV ) ~ 160 MeV |39] , although we 
have considerably increased the uncertainty. We evolve the transversal constant 
in the renormalization scale as 



fK*,±(p) = fK*,±{po) 



a 



4/23 



For the S-meson, only two "moments" are needed, 



duj- 



doo- 



ms $R_(U>) 



m B oj — q z 



16 



(104) 

(105) 
(106) 
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The wave functions $b,±(u>) are modelled such that they fulfill certain relations 
derived from the equations of motion 23 71 . An important property is that 
A^ 1 _(q' 2 ) diverges logarithmically as q 2 — > 0, 



oo, 



(107) 



However, this IR divergence only afflicts the longitudinal (helicity-zero) part of 
the decay amplitude, for which it indicates an enhanced long-distance sensitivity 
at low q 2 . The model employed in Ref. 23 71 leads to the following expressions 
for the 5-meson moments: 



iM), 



-q 2 / (m B Uo)) 



0J 



(-Ei(g 2 / \m B u )) + ire) 



(108) 
(109) 



where u = 0.5(1) GeV. For f B , we use an average of the current LQCD values 
with a generous error bar covering all of them f B = 190(20) MeV [120] . In Table|8 
we collect the numerical values for the hadronic parameters used in this work. 



Table 8: Summary of the K* and B parameters used in this work. 



f K * 


f K *,± (1 GeV) 


a l,a 


Q>2,a 


fB 


^+ 


220(5) MeV 


170(20) MeV 


0.2(2) 


0.1(3) 


190(20) MeV 


2.0(5) GeV" 1 



A. 4 Form factors 

V denotes a vector meson, A its helicity, e\ the corresponding polarization vector, 
etc. B is one of B + , B°, B s . q = p — k. 



(V(k,\)\q lfl b\B(p)) 
(V(k,\)\q W b\B{p)) 



q»(V(k,X)\qa, u b\B(p)) 
q»(V(k,\)\qa, ul5 b\B(p)) 



e e* v n p k u 



i{el-q)%2m v A,{q 2 ) 
+i(m B + m v ) \ e*£ - 



V(q 2 ), 



-a ■ ?y 



;no) 



Mq 2 



{m B -m v Y- A 2 {q\\\) 



rriB+mv 

2ze / ^ pCT e*V^T 1 (g 2 ), (112) 
{e^(m 2 B -m 2 v )-(e*.q)(2p-q),)T 2 (q 2 ) 

+ (^ " ?) U ~ (2p - q),) T 3 (q 2 ). (113) 

V " L B " L V J 
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The matrix element of the scalar density qb vanishes by parity Furthermore, by 
PCAC, 

i{V(k,X)\q l5 b\B(p)) = -^^(e*. q )A (q 2 ), (114) 



(V(k,X)\qP L>R b\B(p)) = ^^_( e *. g L4 (g 2 ). (115) 

m b + m q 

Our sign convention for the Levi-Civita symbol is that of Bjorken and Drell, 
^0123 = +1- Our conventions agree with those of [39]. Some authors, in particular 
[23] , agree on the physical meaning of the form factor symbols (in particular their 
signs), but differ in phase conventions for the hadron states, such that the form 
factor decompositions look different. 



B Effective Lagrangian with vector resonances 

We use the anti-symmetric representation of the vector-meson fields. We explain 
here the particular conventions used in the model and we show some of the re- 
sults for its contribution to the amplitude. Further details on the anti-symmetric 



description of the spin-1 fields can be found in Ref. 97 whereas a discussion 



on the equivalence of this representation with any other one consistent with the 



symmetries and asymptotic behavior of QCD, can be found in Ref. 98 . We de- 
scribe the Si7(3)-multiplets of resonances using flavor-matrices of anti-symmetric 
tensor fields (e.g. V^ v for the nonet of vectors). The interaction Lagrangian for 



resonances 97 98 , at leading order of the chiral expansion in the light sector, is 



4 2) = ^(V^F+) + i^(V^u,u v ). (116) 



In this Lagrangian, u v is the vielbein introducing the pseudoscalar mesons (tt, 
K, rjs) chirally coupled to the resonances fields and F+ u is a generalized external- 
vector (electromagnetic) field strength, that in our case is F+ y ~ 2eQF llu , with Q 
the charge operator in S77(3),F-space. The respective pieces come accompanied 
by the coupling constants fy and gy, which are defined in the chiral limit (the 
same as the mass term in the kinetic Lagrangian). The fields are normalized such 
that 

(0\V^\V(p)) =im v 1 (p"e"(p) -pV(p)) . (117) 
so the vector-resonance fields have the same convention as the one employed in 



the definition of the form factors (see Sec. A.4). Finally, V^ u includes the u and 
</> mesons in ideal mixing, 

oj = -j= {uu + dd) , <p = ss. (118) 
v2 
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In the present case, we need only to consider the LHS of the Lagrangian 



in Eq. (116), with the effective coupling f v receiving small SU(3)p corrections. 
We fix fv for each resonance using the measured V — > e + e~ decay ratios. An 
explicit calculation shows that the decay amplitude of an on-shell resonance V 
with polarization A into an electron (positron) with momentum g_ (g + ) and 
polarization s (s') is 

M v ^ e+e - = i g2/ ^ (iT - q -X] a)7„«(ff+, (119) 
my \ m v J 

where Qy depends on the charge of the constituent quarks in V, Q po = l/v2, 
Q w = l/3y/2 and = -1/3, giving 



^a 2 f v Qy 

V^e+e 7, • 

3my 



(120) 



In Table [9] we list the effective parameters necessary for the description of the 
light-resonance propagation and decay, namely the masses, the decay width, the 
branching fraction for the V — >■ e + e~ decay and the values of fy. For the latter, we 
obtain the errors summing up in quadrature those of the parameters propagated 



linearly through Eq. (120). We have also implemented the S , [/(3)i?-breaking 



explicitly in the meson masses by using the measured Breit-Wigner values 117 . 



Table 9: Hadronic parameters of the light and neutral meson resonances. 



V 


m v [MeV 


Y v [MeV] 


By^ e+e - x 10" 5 


fv [MeV] 


Po 


775.49(34) 


149.1(8) 


4.72(5) 


221(1) 




782.65(12) 


8.49(8) 


7.28(14) 


198(2) 





1019.455(20) 


4.26(4) 


29.54(30) 


228(2) 



On the other hand, the B — > VK* decay amplitude can be written as 
and the contribution of the resonances to the decay is then 



M 



B^rV {-+(.+ £~)K* 



4ira em f v Q^ 



m v (q 2 — fn\ + im v T\ 



:i21) 



:i22i 



where we ignore the off-shell dependence of the real and imaginary parts of the 
pole position and of the B — > VK* amplitude, all of which we expect to have 



in terms of the three helicity amplitudes H v ' oc ey(0, ±)e^*(0, ±)A4 
In principle, 5 different independent observables of the non-leptonic decay are 
required to extract, up to a global phase, the real and complex parts of these 



a negligible effect. The decay amplitude in Eq. (j 12 1 [) is usually paramet erized 



107 
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Table 10: Results for the B — > VK* helicity amplitudes and on the relevant 
observables, which are compared to experimental data 



121 





Po 


u 


<P 




Theo. Expt. 


Theo. Expt. 


Theo. Expt. 


B x 10 6 


2.7(1.2) 3.9(0.8) 


3.3(1.4) 2.0(5) 


9.4(3.8) 9.8(0.7) 


F L 


°A 




31(20) 40(14) 


53(20) 70(13) 


59(19) 48(3) 


h 


o" 




160(40) 


120(31) 


139(20) 136(8) 




-0.11(11) -0.06(9) 


14(15) 0.45(25) 


0.01(5) 



amplitudes. One takes the decay rate plus two out of three polarization fractions, 
h, f±, f\\, and two phases, 



T = {\A \ 2 + \A,_\ 2 + \A^\\ 2 ) , (123) 



± ,H = argA^, (125) 

^0 

where we have used the transversity basis for the amplitudes^ Not all the 
amplitudes are equally important and in fact, in naive factorization, there exists 
a hierarchy among them, 

H° V :H V :H+ = 1: A : ( AY , (126 ) 
m b \m b J 

which results from the V — A nature of the weak interactions and the approximate 
chiral-symmetry of the strong interactions at high energies. However, Hy does 
not factorize in QCDF, although nonfactorizable corrections have been estimated 
(including modelling of end-point- divergent convolutions) and found to be large, 
such that Hy can be comparable to H v . No factorization formula for Hy is 
known, but for B° — > <j)K*° it has been extracted from experiment and found to 
be small, consistent with zero. This reduces to 3 the number of independent ob- 
servables needed for the description of these non-leptonic decays. We summarize 
in the following the results for the two relevant amplitudes H v and Hy in QCD 



factorization 107 and we choose a normalization of the different pieces involved 
that makes the comparison with the other contributions to the semileptonic decay 
more transparent. We also discuss the corresponding predictions for the decay 
rate, F L , 0n and A CP compared to experimental data. 



3 We use a different sign convention for the definition of the transversal amplitude A± as 



compared with Rcf. 107 . 
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The different amplitudes can be written as a sum of products of CKM matrix 
elements, factorizable coefficients containing form factors and decay constants and 
flavor amplitudes containing Wilson coefficients and perturbative corrections in 



the heavy quark limit. In case of the decays under study, we have 107 



* p=c,u * ^ \ 



h b^<pk*° = m B A R »^{a\ + a c 3 ), (127) 



p=c,u 

A 

'A t 

where Aj = VibV*. and we have omitted the polarization index. The valu es fo r 
the amplitudes af\ /5- and d 4 that we use are listed in Table 2 of Ref. 
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An exception is the value for 6i\ that is fixed with the polarization data on the 
B -»■ <f>K* decay (Eq. (39) in Ref. flQ7] ). 



The factorizable coefficients are 



A° VlV2 = ^4*^(0), A± iV2 = ^^-Ff^(O), (128) 

JK* m B JK* 

with ^"^(O) and F±~^ Vl (0) helicity vector form factors in the definition em- 
ployed in Ref. |107 . For the B — > K* form factors we use the values obtained in 



Sec. 3.1[ whereas for the other transitions we take the values given in that refer- 



ence. In Table [10} we show the numerical values for the helicity amplitudes for 
the decays of interest and the respective values for different observables compared 
with experimental data. The errors are obtained adding linearly those from the 
input quantities using Monte Carlo techniques. 

The contribution of the light resonances to the semileptonic decay amplitudes, 
in terms of those of the B — > VK* decay, can be expressed as 

r 0,± a em Gf A t 8lT Qy j K* fv ( m B \ TT 0,± 



= „ 2 2 — " ^ Hp*. (129) 

sl ' v y/2 {q 2 ~m 2 v + im v T v ) \m v J v 

Notice that the amplitude Hy is completely suppressed in the SM so that the 
amplitudes Hy(B — > K*^*) in the semileptonic decays are free from this hadronic 



pollution. On the other hand, and despite of the expected hierarchy in Eq. (126), 
the results in Table 10 show that numerically \Hy \ ~ \H V \, leading to polarization 



fractions that are consistent with the experimental values. 



Finally, and as explained in Sec. 3.2.3, we make the light-quark contribution 



calculated in the VMD model to converge to the one obtained in QCDF around 
q 2 ~ 2 GeV 2 , 

~had,lq = ^had, lq, VMD ( j _ /(g2)) + -had, lq, QCDF /(?2)> (13Q) 
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by using the function 

/fa 2 ) = 1 + exp[ - atf - 2)1 - < 131 > 

Our numerical results are insensitive to the specific choice of a > 0, that we take 
to be a = 12. 
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